Astronomy &: Astrophysics manuscript no. belloche'020712 


February 1, 2008 


(DOI: will be inserted by hand later) 





Molecular line study of the very young protostar 
I RAM 04191 in Taurus: Infall, rotation, and outflow 

A. Belloche^, P. Andre\ D. Despois^, and S. Blinder^'^ 



(N 
O 
O 

(N 



> 

oo 
(N 
l> 
O 
(N 
O 

O 



^ Service d'Astrophysique, CEA/DSM/DAPNIA, C.E. Saclay, F-91191, Gif-sur-Yvette Cedex, France 
^ Observatoire de Bordeaux (INSU/CNRS), B.P. 89, F-33270 Floirac, France 

Division of Nuclear Medicine, Vancouver Hospital and Healtii Sciences Center, Vancouver, B.C., Canada 

Received April 16, 2002; accepted July 12, 2002 

Abstract. We present a detailed millimeter line study of the circumstellar environment of the low-luminosity 
Class protostar IRAM 04191 + 1522 in the Taurus molecular cloud. New line observations demonstrate that 
the ~ 14000 AU radius protostellar envelope is undergoing both extended infall and fast, differential rotation. 
Radiative transfer modeling of multitransition CS and C'^^S maps indicate an infall velocity Dinf ~ 0.15 km s~^ 



at r ~ 1500 AU and 



"inf 



0.1 km s up to r ~ 11000 AU, as well as a rotational angular velocity Q 



at ~ 11000 AU. Two distinct regions, which differ in both their infall and their rotation properties, therefore 
seem to stand out: the inner part of the envelope (r ^ 2000 — 4000 AU) is rapidly collapsing and rotating, while 
the outer part undergoes only moderate infall/contraction and slower rotation. These contrasted features suggest 
that angular momentum is conserved in the collapsing inner region but efficiently dissipated due to magnetic 
braking in the slowly contracting outer region. We propose that the inner envelope is in the process of decoupling 
from the ambient cloud and corresponds to the effective mass reservoir 0.5 Mq) from which the central star 
is being built. Comparison with the rotational properties of other objects in Taurus suggests that IRAM 04191 
is at a pivotal stage between a prestellar regime of constant angular velocity enforced by magnetic braking and 
a dynamical, protostellar regime of nearly conserved angular momentum. The rotation velocity profile we derive 
for the inner IRAM 04191 envelope should thus set some constraints on the distribution of angular momentum 
on the scale of the outer Solar system at the onset of protostar/disk formation. 
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1. Introduction 

1.1. Tiie eDigiudl'iL uiii,eL of pruLu^>Lelldr culldfji^e 



accretion history if the collapse initial conditions are ei- 
ther not singular or not scale-free (e.g. Foster & Chevalier 



1993; Henriksen, Andre fc Bontcmps 1997; Basu 1997 



Despite recent progress, the initial conditions of star for 
mation and the first phases of protostellar collapse remain 
known (e.g. Myers 199?)[ [Andre, Ward-Thompson 



poorly 



[Ciolek fc Konigl 1998| ; [Hennebelle et al. 2002| ). Starting 
from realistic, finite-sized prestellar cores with p « cte in 
their central region (c.f. Ward-Thompson, Motte fc Andre 



fc Rar Lnv 2nfin| fofl^^^^d^^^^Tn the^^TTT^ 1999| ; pacmann et al. 200q ; [Alves, Lada fc Lada 2001| ) 



isolated, low-mass star formation (e.g. Shu, Adams fc 



Lizano 1987), the initial conditions correspond to essen- 



tially static singular isothermal spheres (SISs, which have 
p = (as^/27rG)r^^), assumed to be in slow, solid-body ro- 

TSC84). This leads 



tation (Terebey, Shu fc Cassen 1984 



to a strictly constant mass accretion rate, Ma, 



(where is the isothermal sound speed), and to a growth 
of the centrifugal disk as Rdisk oc (cf. |TSC84 ) during 

the prntnstpllar a.rrrptinn phasp (/ -> 0) Othnr thpnrpli- 



these models yield supersonic inward velocities close to 
the center prior to point mass formation (i.e. at t < 0) 
and result in denser, nonequilibrium density distributions 
with strong differential rotation at the onset of the main 
accretion phase, i.e., at t = 0. In these models, the accre- 
tion rate Mace is initially significantly larger than in the 
Shu model, then quickly converges toward the standard 
~ aJ^/G value, and finally declines much bel ow Qs^/G 



because of the finite reservoir of mas s (see, e.g., [Foster fc 



cal mo Hels pyist, howPver, that prpdict a timp-dpppudpnt Chevalier 1993[ ; [Henriksen et al. 1997[ ). Conservation of an- 
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gular momentum during dynamical collapse at t < pro- 
duces a differential rotation profile at t = (e.g. cx r^^ 
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in the magnetically-controlled model of Basu 1998). This 
rotation profile in turn implies a more rapid growth of 
Rdisk initially (i.e., Rdisk oc i at small i > in the Basu 



model) than in the |rSC84| model. 

Getting at a better, more quantitative knowledge 
of protostellar collapse is crucial, e.g., to gain insight 
into the origin of stellar masses and disk formation. 
Observationally, there are two complementary approaches 
to estimating the initial conditions of protostar formation. 
The first approach consists in studying the structure and 



kinematics of "prestellar cores" such as L1544 (e.g. Ward- 



Thomf son et al. 1999; Tafalla et al. 1998) , representative 
of times t ^Q. The second approach, adopted here, is the 
detailed study of Class accreting pro tostars observed at 



t > 0, such as IRAM 04191 (see § |lj), which should still 



retain detailed memory of their initial conditions. 

1.2. IRAM 04191: A very young Class protostar in 
Taurus 

The massive {Mtot ^ 1.5 Mq) dense core/envelope of the 
Class object, IRAM 04191 + 1522 (hereafter IRAM 
04191), was originally discovered in the millimeter dust 
continuum with the IRAM 30m telescope in the south- 
ern part of the Taurus molecular cloud ( [Andre, Motte 
fc Bacmann 1999| - hereafter AMB99). Follow-up obser- 



vations revealed a highly coUimated CO bipolar outflow 
(see Fig. |l|), a weak 3.6 cm VLA radio continuum source 
located at its center of symmetry, and spectroscopic evi- 
dence of spatially extended infall motions in the bulk of 
the env elope. These are typical attributes of a Cla s s pr o- 
tostar ([Andre, Ward-Thompson fc Barsony 1993| , |2000|) . 



The very high envelope mass to luminosity ratio of 
IRAM 04191 (M<f > 3 Mq/Lq) and its posi- 

tion in the Lboi — ^boi evolutionary diagram (Lboi ~ 0.15 
L0 and Tboi ~ 18 K) suggest an age t ~ 1— 3x 10* ^ yr since 
the beginning of the accretion phase (see AMB99| ) . This is 
significantly younger t han all of the IRAS candidate pro 
tostars of Taurus (e.g. Kenyon, Calvet fc Hartmann 1993[) 



including L1527 which has t ^ 10 yr (e.g. Ohashi et al 



1997a ). IRAM 04191 thus appears to be the youngest ac- 
creting protostar known so far in Taurus, although the 
collapsing protostellar condensa tion MC 27 discovered by 
Onishi, Mizuno fc Fukui (1999| ) may be in a comparable 



evolutionary state. 

As IRAM 04191 is particularly young, nearby (d = 
140 pc), and relatively isolated, the study of its velocity 
structure based on molecular line observations provides a 
unique opportunity to set constraints on collapse models. 
This is especially true since the viewing angle is favor- 
able. The C0(2-l) outflow map of |AMB99| (see Fig. ^ 
shows well separated outflow lobes with almost no overlap 
between blue-shifted and red-shifted emission, indicating 
that the flow lies out of the plane of the sky at an inter- 
mediate inclination angle (e.g. Cabrit fc Bertout 1990| ). 
In addition, [AMB99| estimate an aspect ratio of ~ 0.65 
for the circumstellar dust/N2H+ envelope, whose major 
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Fig. 1. N2lI+(l-0) integrated intensity map of the 
IRAM 04191 protostellar envelope overlaid on the C0(2- 
1) outflow map of [AMB9E . Both maps were taken with 
the IRAM 30m telescope. The N2lI+(l-0) emission is in- 
tegrated over the whole seven-component multiplet from 
5.7 to 22.6 km s~^ and the contours go from 0.5 to 4.5 
K km s^i by 0.5 K km s'^. The blueshifted C0(2-l) emis- 
sion {dotted contours) is integrated over the 0-5 km s~^ ve- 
locity range, and the redshifted C0(2-l) emission {dashed 



contours) is integrated over the 8-13 km 



velocity 



range; the C0(2-l) contours go from 5 to 21 K km s ^ 
by 2 K km s^^. 

The insert in the bottom-left corner shows a 227 GHz 
dust continuum map of the inner envelope o btain ed at 
the IRAM Plateau de Bure interferometer (see § 3T), with 



positive (solid) contours from -1-1 to -1-6 by 1 mJy/1.9"- 
beam; the dotted contour is negative at —1 mJy/1.9"- 
beam). 

In both maps, the central star symbol at (0,0) marks the 
position of IRAM 04191 as determined by a 2D Gaussian 
fit to the PdBI 227 GHz continuum image. The other 
star symbol indicates the position of the Class I source 
IRAS 04191. 



axis is perpendicular to the outflow axis (cf. Fig. [l|). Both 
characteristics are consistent with an inclination angle of 
the outflow axis to the line of sight of z ^ 50°. 

Here, we present and discuss the results of a com- 
prehensive set of molecular line observations toward 
IRAM 04191. The layout of the paper is as follows. Sect. | 
summarizes observational details. Sect. ^ interprets the 
observations in terms of infall, rotation, and outflow mo- 
tions in the envelope. We then model the observed spectra 
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using radiative transfer simulations computed in ID spher- 
ical geometry with radial infall motions (Sect. |^) and in 
2D axial geometry with both infall and rotational motions 
(Sect. |). Sect. I compares the derived constraints on the 
velocity structure of the IRAM 04191 envelope with the 
predictions of collapse models. Our conclusions are sum- 
marized in Sect. 0. 

2. Observations 

We carried out millimeter line observations with the 
IRAM 30m telescope at Pico Veleta, Spain, during 7 
nights in July and August 1999, in the following molecular 
transitions: N2H+(l-0), CS(2-1), C34S(2-1), HCO+(1-0), 
H13CO+(1-0), C3H2(2i2-loi) at 3 mm, CS(3-2), C^'^S{3- 
2), H2CO(2i2-lii) at 2 mm, and CS(5-4), H2CO(3i2- 
2ii), 01^0(2-1), DCO+(3-2) at 1.2 mm. Our adopted set 
of line frequencies is given in Table |l|. The half-power 
beamwidths were ^ 26 ", ^ 17 " and ~ 10 " at 3 mm, 
2 mm and 1.2 mm, respectively. We used four SIS het- 
erodyne receivers simultaneously and an autocorrelation 
spectrcjmeter as backend, with a spectral resolution ot 



Table 1. Adopted line rest frequencies and telescope effi- 
ciencies. 



Line 


Frequency'^' 




Ref. 








(MHz) 


'km s""'" 


) (3) 


(%) 


(%) 


C3H2(2l2-loi) 


85338.905(6) 


0.02 


(1) 


92 


73 


N2H+(ioi-oi2) 


93176.258(7) 


0.02 


(2) 


92 


73 


C^^S(2-1) 


96412.952(1) 


0.003 


(3) 


92 


73 


CS(2-1) 


97980.953(1) 


0.003 


(3) 


90 


73 


C^*O(l-0) 


109782.175(2) 


0.005 


(4) 


92 


73 


C^^S(3-2) 


144617.101(1) 


0.002 


(3) 


87 


65 










90 


54 


CS(3-2) 


146969.026(1) 


0.002 


(3) 


90 


54 


Ci**0(2-1) 


219560.3541(15) 


0.002 


(5) 


86 


42 


CS(5-4) 


244935.555(1) 


0.001 


(3) 


84 


42 



Notes: "-^^ The frequency uncertainty in units of the last sig- 
nificant digit is given in parentheses. 
Frequency uncertainty converted in units of velocity. 
References fo r rest frequencies: (1) Laboratory 



(2) 
(3) 



measure- 



Vrtilek, 



ment from 
Observational result 
(3) Laboratory measurement 



Gottlieb fc Thaddeus (1987|); (2 ) 
from [Lee, Myers fc Tafalla (200l"|) 



20 kHz at 3 mm and 40 kHz at 2 mm and 1.2 mm. i'hc cor- 
responding velocity resolutions ranged from 0.05 to 0.08 
km s~^ per channel. The observations were done in sin- 
gle sideband mode, with sideband rejections of 0.01, 0.1 
and 0.05 at 3 mm, 2 mm and 1.2 mm, respectively. The 
forward and beam efficiencies of the telescope used to con- 
vert antenna temperatures Tf into main beam tempera- 
tures Tmb are listed in Table 0. The system temperatures 
(TJ scale) ranged from ~ 110 K to 150 K at 3 mm, 
~ 280 K to - 410 K at 2 mm, and - 300 K to - 550 
K at 1.2 mm. The telescope pointing was checked every 
- 2 hours on Saturn, 0528+134, and/or 0420-014 and 
found to be accurate to ~ 3 " (rms). The telescope fo- 
cus was optimized on Saturn every ~ 3 hours. Position 
switching observations were done with a reference posi- 
tion located at either {Aa,AS) = (1200 ",-1200 ") or 
(80 ",—80 ") relative to the (0,0) position (envelope center 
as measured in the 1.3mm continuum). E xtensive mapping 
was p erformed in the "on-the-fly" mode ( Ungerechts et al 
^000| ). A few additional Ci**0(2-1), CS(3-2), 

C34S(3-2), CS(5-4) and Hi3CO+(3-2) position-switch ob- 
servations, performed in September 1997, November 2000 
and October 2001, will also be used here. All of these 
single-dish data were reduced with the CLASS software 
package (Buisson et al. 2002). 

In addition, we also observed IRAM 04191 with the 5- 
antenna IRAM Plateau de Bure interferometer (PdBI) in 
its Bl, C2, D configurations between January and April 
1999. The two receivers of each antenna were tuned to 
the CS(2-1) and H2CO(3i2-2ii) lines, with spectral reso- 
lutions of 40 kHz and 80 kHz, corresponding to velocity 
resolutions of 0.12 km s^^ and 0.10 km s""'^, respectively. 
The four remaining windows of the PdBI correlator were 
used to record the continuum emission with a total band- 
width of 300 MHz at both 98 GHz (A ~ 3 mm) and 
227 GHz (A ~ 1.3 mm). The (naturally- weighted) syn- 



Thaddeus (2002). (4) 



from 
Laboratory 



Gottlieb, Myers & 



measurement from 
pr ivate communication ) . (5) Laboratory 



Klapper (2001, 
measurement from 
Forward and beam efficiencies of the IRAM 30m telescope 



Klapper et al. (2001) 



(Wild (1999) and tittp://www.iram.es/) 



thesized half-power beamwidths were 4.5 " x 4.4 " (630 
AU X 620 AU) at 98 GHz and 1.9 " x 1.8 " (270 AU x 
250 AU) at 227 GHz, and the (FWHP) primary beams 
^ 50 " and ^ 25 ", respectively. The correlator band- 
pass was calibrated on the strong source 3C273. Several 
nearby phase calibrators were observed to determine the 
time-dependent complex antenna gains. The absolute cal- 
ibration uncertainty was estimated to be ~ 15%. The data 



were calibrate d and imaged using the CLIC (Lucas 1999 ) 
and Mapping (Guilloteau, Lucas & Bouyoucef 2002) pack- 
ages in the IRAM software. 

3. Results and qualitative interpretation : 
Evidence for rotation and infall 

3.1. Weak 1.3 mm continuum detection at PdBI 

The PdBI continumn image at 227 GHz (see Fig. ^ 
reveals weak, point-like emission centered at Q!(2ooo) = 
04'^21'»56"91, (5(2000) = 15°29'46.1", a position which 
should be accurate to better than 0.5". (The PdBI po- 
sition is offset by ^ 5 " from the N2H"'"(l-0) centroid ob- 
served at the 30m telescope, which is only marginally sig- 
nificant given the 3 " single-dish pointing accuracy.) 

1 q" 

We measure a peak 227 GHz flux density of S'p^ak — 
6.1 ± 0.4 mJy/1.9"-beam at the central position. The 
227 GHz emission is slightly resolved and a Gaussian fit 
performed in the uv-plane yields a deconvolved FWHM 
size of (1.1" ± 0.4") x (0.6" ± 0.3") with a position an- 
gle P.A. = 84 °. At 98 GHz, the rms noise level is 0.14 
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inJy/4.5"-beani and we do not detect any emission above 
0.6 mJy/4.5"-beam (- 4ct)[|. 

For comparison, using the IRAM 30m telescope 
equipped with the MP IfR bolometer array (MAMBO), 
Motte & Andre (2001) measured a peak flux density of 

'5'poak = 110 ± 7 mJy/ll"-beam at 1.3 mm {veff 
240 GHz - e.g. Broguicrc ct al. 2002) and a radial in- 



tensity profile of the form I {6) cx 9 0-6±o.i jj-^ ^j-^g range 
of angular radii 6* = 11 " to 6* = 100 ". This extended 
1.3 mm continuum source is clearly the dust counterpart of 
the circumstellar gas envelope/core observed in N2H+ (cf. 
Fig. |l|). Assuming the same radial intensity profile holds 
at smaller angular radii, one expects the peak flux den- 
sity of the envelope to scale as S'pcak(^b) oc i with 
beamsize Ob- If we adopt a dust opacity index (3 — 1.5 to 
account for the slight difference in observing frequency be- 
tween the 30m and PdBI measurements, this flux-density 
scaling predicts S^ f^^. ~ 7.7 mJy/1.9"-beam, which is only 
25% larger than the PdBI peak flux density quoted above. 
Given the relative calibration uncertainties, this compar- 
ison suggests that the weak 227 GHz emission detected 
at PdBI arises from the inner part of the envelope seen 
at the 30m telescope rather than from an accretion disk 
surrounding the central protostellar object. Furthermore, 
a single power law I{9) oc Q-0-5±0-2 appears to charac- 
terize the radial intensity profile of the envelope over the 
whole range of angular radii from ~ 1" to ^ 100". Such 
an intensity profile implies either a relatively fiat p cx r"^ 
density profile if central heating with T^ust r~'^'^ ap- 
plies in the inner part of the envelope (see Fig. ^ below), 
or a p cx r^^-^ density profile if the dust temperature is 
uniform. 

The weak emission detected at PdBI sets strong con- 
straints on the mass and size of any central accretion disk. 
Assuming optically thin emission, a dust mass opacity of 



0.02 cm^ g ^ typical of circumstellar disks (cf. Bcckwith 
et al 

1 



.990) and a mean dust temperature of 20 K, the 
PdBI peak flux density yields Mjisk < 1 x 10"^ Mq. 



< 1 X 10' 

Alternatively, the assumption of optically thick dust emis- 
sion at 20 K implies i?disk < 10 AU. 

3.2. Fast, differential rotation 



All of the centroid velocity (first-order moment) maps 
taken at the 30m telescope in small optical depth lines, 
such as C3H2(2i2-loi), N2H+(101-012), Hi3cO+(l-0), 
and C'^^S(2-1), show a clear velocity gradient across the 
envelope (see, e.g., Fig. ||). The (south)-east part is red- 
shifted with respect to the source systemic velocity, while 
the (north)-west part is blueshifted. We have applied the 
sector method described by Arquilla & Goldsmith (1986) 



to measure the direction of this velocity gradient in the 
C3H2(2i2-loi) rnap using Gaussian fits to the spectra. We 



^ We measure a peak 98 GHz flux density of 0.6 mJy/4.5"- 
beam close to (~ 4" north-east of) the central 227 GHz posi- 
tion but the cleaned image contains negative contours down to 
-0.8 mJy/4.5"-beam. 
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Fig. 2. Map of the peak LSR velocity {filled squares) in 
the IRAM 04191 envelope as derived from Gaussian fits 
to the C3H2(2i2-loi) spectra observed with the 30m tele- 
scope. The map of integrated C3H2(2i2-loi) intensity be- 
tween 5.7 and 7.7 km s~^ is overlaid as contours (ranging 
from 0.4 to 1.0 by 0.2 K km s^^). The black arrow shows 
the best- fit direction (P. A. 114° ± 8°) of the mean ve- 
locity gradient {^-^ 3.4 km s~^ pc~^) across the envelope; 
the grey-shaded sectors display the la error on this orien- 
tation. The size of each square is proportional to the LSR 
velocity at its position. The open squares in the upper- 
left corner give the LSR velocities measured at positions 
that were ignored in fitting the velocity gradient because 
of likely contamination by the Class I source IRAS 04191 
and/or the redshifted lobe of the IRAM 04191 outfiow. 
The two star symbols mark the positions of IRAM 04191 
[at (0,0)] and IRAS 04191. 



obtain a position angle P. A. ~ 114° ±8° (cf. Fig. |). The 
velocity gradient thus lies along the major axis of the elon- 
gated dust/N2H+ core (P. A. ^ 120°), i.e., perpendicular 
to the outfiow axis. Coupled to the high degree of symme- 
try of the position- velocity diagrams shown in Fig. ^ with 
respect to envelope center, this strongly suggests that the 
envelope is rotating about an axis coinciding with the out- 
fiow axis (cf . Fig. nl) . Turbulent motions would produce a 



more random velocity field (cf. Burkcrt & Bodenheimei 
2OOOD and are weak here anyway (see § 3.4 and Fig. |l2|b 



below). 

From now on, we adopt a position angle P.A. = 28° for 
the projection of the rotation/outflow axis onto the plane 
of the sky. To minimize contamination by the outfiow, 
we analyze the velocity structure of the envelope along 
the axis perpendicular to the outflow (P.A. = 28°) and 
going through the center. Along this axis, the magnitude 
of the velocity gradient estimated from linear fits to the 



^ Before fitting the direction of the velocity gradient, we have 
masked the redshifted emission seen in the north-east part of 
the C3H2(2i2-loi) map shown in Fig. |, probably due to the 
redshifted lobe of t he outflow or to the nearby Class I source 
IRAS 04191-^1523 ([Tamura et al. 199l|). 
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Fig. 3. Position-velocity diagrams in small optical depth 
lines along the direction perpendicular to the outflow axis, 
(a) C3H2(2i2-loi); (b) Hi3CO+(l-0); (c) C34S(2-1); (d) 
N2H+ (101-012). The dots with error bars represent the 
peak velocities as determined by Gaussian fits. The base 
contour and contour step have the same value: 0.3 K (in 
T;^) in (a) and (b); 0.2 K in (c) and (d). Note how the dots 
mark a similar S shape in all four diagrams, suggestive of 
differential rotation in the envelope (see text). 



C3H2(2i2-loi), Hi3CO+(l-0), m{2-l), and N2H+(101- 
012) position- velocity diagrams is ~ 6 km s^^ pc^^, ^ 
7 km s~^ pc~^, ~ 6 km s~^ pc~^, and 10 km s~^ pc^^, 
respectively, increasing from north-west to south-east over 
40 " (see Fig. |^). The mean velocity gradient is 7 ± 2 
km s~^ pc~^. Given the estimated viewing angle i — 50° 
of the flattened envelope (cf. § ^^), this implies a mean 
rotational angular velocity of 17 = 9 ± 3 km s^^ pc^^ in 
the inner r ^ 2800 AU radius region^. 

The C3H2(2i2-loi), Hi3CO+(l-0), and N2H+(101- 
012) position-velocity diagrams shown in Fig. ^ may be 
slightly contaminated by a secondary, redshifted compo- 
nent toward t he east. We have used th e GAUSSCLUMPS 
algorit hm o f 3tutzki fc Giistcn (1990| ) (see also Kramer 
et al. ^99?) to try and subtract this secondary compo- 
nent. The algorithm finds a Gaussian component at ap- 
proximately (-1-50 ", +30 ") in C3H2(2i2-loi) and (-1-30", 
+2Q") in H"CO+(1-0). The removal of this secondary 
component reduces the magnitude of the velocity gradient 
measured in C3H2(2i2-loi) and Hi3CO+(l-0) by - 30%. 
GAUSSCLUMPS fails to identify any secondary compo- 



^ With such an angular velocity, the gas at r ~ 2800 AU 
would make a complete turn in ~ 0.7 Myr, which is compara- 
ble to the typical lifeti me of prestellar cores with centra l densi- 
10^ cm"^ (e.g. llessop & Ward- Thompson 2000|). Thus, 
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although the observed velocity gradient is probably not indica- 
tive of a well-developed circular motion, it may have already 
induced significant rotational distortion in the envelope. 
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Fig. 4. Normalized C3H2(2i2-loi) and C-'^'*S(2-1) spectra 
observed every half beamwidth along the direction perpen- 
dicular to the outflow axis. Each spectrum has been di- 
vided by its peak temperature for easier inter-comparison. 
In both panels, the vertical dotted line marks the peak 
velocity of the central spectrum {dash- dotted curve at the 
bottom) measured by a Gaussian fit. The double-sided 
arrows and associated values indicate the magnitude of 
the velocity shifts (determined by Gaussian fits) between 
spectra observed at symmetric positions (±13", ±26", 
±39", or ±52") about the center. The channel resolution 
is shown in the upper right corner of each panel. While 
solid-body rotation would produce a velocity shift increas- 
ing linearly with radius, the opposite trend is observed be- 
yond one beamwidth (26") from the center: the velocity 
shift decreases for larger offsets about the central position. 
This strongly suggests the presence of differential rotation 
in the envelope. 

nent near the same location in the N2H+ (101-012) data. 
The velocity gradient thus seems to be significantly larger 
in N2H+ (101-012) than in the other three lines (see be- 
low). 

Remarkably, all of the position-velocity diagrams 
shown here present a S shape, clearly seen on the peak 
velocity curves derived from Gaussian fits (see Fig. ||). 
On either side of the source, the velocity shift with re- 
spect to the systemic velocity increases in absolute value 
up to an angular radius ^ 25 ± 10" (i.e., one full 
beamwidth of the 30m telescope at 3mm, corresponding to 
3500±1400 AU), and then decreases for 61 > 6',„. This 
is further illustrated in Fig. ^ which shows the C3H2(2i2- 
loi) and C^'*S(2-1) spectra observed every half beamwidth 
(Nyquist sampling) at symmetric positions with respect to 
source center along the direction perpendicular to the out- 
flow axis. It can be seen that the velocity shift between 
symmetric positions decreases from 0.17 km s~^ (i.e., 2.5 
channels) at ±26" to 0.09 km s^^ (1.3 ch.) at ±52" in 
C3H2(2i2-loi). In C34S(2-1), the velocity shift varies from 
0.16 km s-^ (2.6 ch.) at ±26" to 0.14 km s'^ (2.3 ch.) 
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at ±39". These values translate into a decrease of the 
C3H2(2i2-loi) velocity gradient from 4.8 km s"^ pc~^ to 
1.3 km s~^ pc~^ between 26" and 52", and to a decrease 
of the C'^^S(2-1) velocity gradient from 4.5 km s~^ pc~^ 
to 2.6 km s^^ pc^^ between 26" and 39". The decrease of 
the velocity shift with radius beyond ±26" is in marked 
contrast with the linear increase expected in the case of 
solid-body rotation. We thus conclude that there is strong 
evidence for differential rotation in the envelope beyond 
~ 3500 AU. 

Correcting for inclination and taking all four lines of Fig. ^ 
into account, we estimate that the rotational angular ve- 
locity decreases by a factor of ~ 5 from = 9 ± 3 
km s-i pc-i (- 3 X 10~i3 rad s'^) at r = 2800 AU (20") 
to n = 1.9 ± 0.2 km s-i pc'^ (- 6 x IQ-^^ rad s'^) 
at r = 7000 AU (50"). At a radius r - 11000 AU, 
our C^^O(l-O) and C^*0(2-l) observations suggest an 
even smaller angular velocity, ^ 0.5 — 1 km s~^ pc~^ 
(~ 1.5 — 3 X 10"^"* rad s~^). On scales smaller than the 
beam, the intrinsic rotation velocity pattern is uncertain 
due to insufficient spatial resolution (the beam HWHM 
angular radius is ^ 13 " at 3 mm, corresponding to a 
physical radius of ^ 1800 AU). However, two indirect ar- 
guments suggest that the differential rotation pattern ob- 
served here between ~ 3500 AU and ~ 11000 AU con- 
tinues down to smaller (~ 1000 — 2000 AU) scales. First, 
such a differential rotation pattern, combined with a lower 
level of molecular depletion near envelope center in N2H"'' 
(cf. § 3.4 and § |4.2| below), would explain the higher ve- 
locity gradient (10 km pc~^) measured in N2H+(101- 
012) over ±20" compared to the other three lines shown 
in Fig. ^. Second, using NH3 interferometric observations 
sensitive to ^ 5" — 15" scales, Wootten, Wiseman & 



Fuller J2OOI) have recently reported an even larger gra- 
dient 15 km s^^ pc^^) than our present N2H+ value. 



3.3. Spectroscopic signature of collapse 



As already pointed out by AMB99, the classical spec- 
troscopic signature of infal l motions (cf. pvans 1999| ; 
Myers, Evans fc Ohashi 200C| ) is seen toward IRAM 04191. 
Optically thick lines such as CS(2-1), CS(3-2), CS(5-4), 
H2CO(2i2-lii) and H2CO(3i2-2ii) are double-peaked and 
skewed to the blue, while low optical depth lines such 
as C34S(2-1) and 0^48(3-2) peak in the dip of the self- 
absorbed lines (see spectra observed at the central posi- 
tion in Fig. I). Blue-skewed CS(2-1) and CS(3-2) spectra 
are observed in an extended region, up to an angular ra- 
dius of at least 40" (5600 AU) from source center (see cut 
taken perpendicular to the outflow axis in Fig. || below). 
Such asymmetric line profiles with a blue peak stronger 
than the red peak are expected in a collapsing envelope 
when the line excitation temperature increases toward the 
center. There is therefore strong evidence for the presence 
of extended inward motions in the IRAM 04191 envelope. 

The blue-to-red peak intensity ratio of the CS lines 
is weaker toward the south-east (i.e. the envelope hemi- 
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Fig. 5. Spectral line profiles observed at the central posi- 
tion of IRAM 04191 in the optically thick CS(2-1), CS(3- 
2), CS(5-4), H2CO(2i2-lii), H2CO(3i2-2ii) lines, and the 
low optical depth 0^48(2-1), C34S(3-2), H"CO+(1-0), 
C3H2(2i2-loi), N2H+(101-012) lines. The dotted line indi- 
cates the source systemic velocity estimated from a seven- 
component hyperfine structure fit to the N2H+(l-0) mul- 
tiplet. 



sphere red-shifted by rotation), while the asymmetry is 
stronger toward the north-west (i.e. the hemisphere blue- 
shifted by rotation) where the red peak is even barely 
visible. This behavior is in qualitative agreement with the 
expected distortion of the infall asymmetry due to rota- 
tion when the rota tion velocit y does not dominate over 
the infall velocity| (|Zhou 1995|) . 

The dips of the optically thick CS(2-1), CS(3-2), and 
CS(5-4) lines have velocities of 6.75 ± 0.03 km s^^, 6.81 ± 
0.04 km s""'^, and 6.86 ± 0.05 km s~^, respectively, and 



are redshifted relative to the source systemic velocity of 
6.66 ± 0.03 km s'^ (see Fig. |). The latter value results 
from a 7-component Gaussian fit to the N2H+(l-0) mul- 
tiplet, using the hyperfine structure (hfs) method of the 
CLASS reduction software, and assuming the relative fre- 
quencies and intensities of the 7 hyperfine components 
determined by [Caselli et al. (1995|), and the N2H+ (101- 
012) rest frequency of |Lee, Myers fc Tafalla (200l[ ) (see 
Table El). The uncertainty on the Lee et al. (2001) fre- 



quency is estimated to be ~ 0.02 km s^ , to which we 
conservatively add (in quadrature) an uncertainty of 0.015 
km s~^ arising from a maximum pointing error of ^ 5", 
given the velocity gradient discussed in § above. This 
gives a final uncertainty of ^ 0.03 km s^^ on the source 
systemic velocity. The dips of the self-absorbed CS(2-1), 
CS(3-2), and CS(5-4) spectra are thus redshifted relative 
to the systemic velocity by 0.09±0.04 km s"\ 0.15±0.05 



When the rotation velocity becomes comparable to the in- 



fall velocity, the asymme t ry may even be reversed (e.g. Walker 
Narayanan fc Boss 199-^ ; [Ward- Thompson fc Buckley 2001| )~ 
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km s~^, and 0.20 ± 0.06 km s~^, respectively. These ab- 
sorption dips are presumably produced by the outer lay- 
ers of the front hemisphere of the envelope. The fact that 
they are redshifted provides a second indication that in- 
ward motions are present in the outer envelope where the 
opacity of the CS lines is of order unity. More quantita- 
tively, the main-beam brightness temperature of the dip 
in the central CS(2-1) spectrum is ~ 0.7 K, which matches 
the peak temperature of the spectra taken at 70 " from 
source center. These spectra are still optically thick since 
we measure a CS(2-1) to C'^'*S(2-1) integrated intensity 
ratio of only ^ 5, i.e., ~ 4 times less than the standard 
CS to S isotopic ratio of ~ 22 ([Wilson fc Rood 1994|). 



Assuming a spherically symmetric envelope, we conclude 
that the absorbing shell producing the dip in the central 
CS(2-1) spectrum has an angular radius larger than ^ 70". 
The observed ^0.1 km s^^ redshift of the dip is thus sug- 
gestive of inward motions 0.1 km extending up to a 
radius of at least 10000 AU. Radiative transfer simulations 
confirm this conclusion (see Sect. ^ below). 

On the other hand, the CS(5-4) emission is much more 
concentrated spatially than the CS(2-1) emission, and con- 
fined to the inner 20 " (FWHM) region (see the non- 
detection at 10" in Fig. || below). The radius of the shell 
producing the absorption dip in the central CS(5-4) spec- 
trum must therefore be smaller than 10 ". The observed 
^ 0.2 km s~^ redshift of the CS(5-4) dip then suggests 
faster inward motions in the central (r < 1400 AU) re- 
gion. 



an optical depth of ~ 0.85 for the isolated N2H+ (101-012) 
component, suggesting neghgible ( ^ 15%) optical depth 
broadening. The level of optical depth broadening should 
thus be more pronounced in C'^'^S and cannot explain the 
difference in linewidth between N2H+ and C'^'^S. 

We propose that this difference in linewidth between 
N2H"'" and C'^'^S results from a combination of higher 
infall/rotation velocities and lower N2H"'" depletion to- 
ward the center. Bergin & Langer (1997) have shown 
that sulphur-bearing molecules such as CS are strongly 
depleted when the density increases, whereas N2H+ re- 
mains in the gas phase, at least up to densities ^ 10^ 
cm^'^. Indeed, we measure a decrease of the C^''S(2- 
1)/N2H"'" (101-012) integrated intensity ratio by a factor 
of 2 from ~ 5000 AU to - 2000 AU. As both hues have 
nearly the same critical density and are approximately op- 
tically thin|^, the decrease of the integrated intensity ratio 
may be interpreted as a decrease of the C'^^S/N2H+ abim- 
dance ratio toward the center. 

Finally, we note that the N2H+(101-012) linewidth 
peaks at the central position (as shown by a linewidth- 
position plot along the direction perpendicular to the on- 
flow axis). N2H"'" may thus be more sensitive to higher 
velocity material produced by, e.g., infall, rotation, or 
outflow near the central protostellar object. As N2H"'" 
is generally underabundant in molecular outflows (e.g. 
Bachiller & Perez- Gutierrez 1997), the central broadening 
of N2H+(l-0) is most likely due to infall and/or rotation 
motions. 



3.4. Linewidths: Evidence for turbulent infall ? 

The C34S(2-1) and C^^S{3-2) spectra displayed in Fig. | 
are slightly asymmetric and skewed to the blue, which sug- 
gests they are marginally optically thic k and showing som e 
infall asymmetry (see, e.g.. Fig. 1 of Myers et al. 1995| ). 
We measure (FWHM) linewidths of 0.38 ± 0.01 km s^^ 
for C34S(2-1) and 0.32 ±0.03 km s'^ for C34S(3-2), which 
are 3.8 and 3.2 times larger than the thermal broaden- 
ing of C'^^S at a kinetic temperature Tk = 10 K, respec- 
tively (see § 4.1 for constraints on Tk). Radiative transfer 
simulations (cf. Sect. ||) indicate opacities of ~ 1.5 and 
- 1 for C34S(2-1) and C34S(3-2), respectively Therefore, 
line saturation effects cannot broaden the C'^''S(2-1) and 
C34S(3-2) spectra by more than - 30% and ~ 20%, re- 
spectively, and the linewidths are primarily nonthermal. 
Motions such as infall, rotation, outflow, or "turbulence" 
along the line of sight are required to explain such nonther- 
mal linewidths. The nonthermal motions do not dominate 
over thermal motions, however, since the C^''S(2-1) and 
C'^''S(3-2) linewidths represent only ~ 80% of the thermal 
velocity dispersion for a mean particle of molecular weight 
H = 2.33. The N2H+ (101-012) line is slightly broader, 
with a FWHM ~ 0.55 km s^^, i.e., 1.2 times broader 
than the (mean particle) thermal velocity dispersion and 
~ 1.5 times larger than the C^'*S linewidths. The hyper- 
fine structure fit to the N2H+ multiplet (see § 3.3) yields 



3.5. The CS line wing emission 

The morphology of the single-dish CS(2-1) maps inte- 
grated over the 4.6-6.1 km s~^ (blue) and 7.1-8.6 km 
(red) velocity ranges strongly suggests that the CS line 
wing emission is dominated by material associated with 
the outflow (see Fig. |a & Fig. ||b). Some CS(2-1) emission 
is detected at the edges of both the red and the blue lobe 
of the CO outflow: the CS(2-1) red wing is relatively weak 
and concentrated at the south-eastern edge and the tip of 
the red CO lobe (Fig. gi), while the CS(2-1) blue wing is 
much stronger and distributed in two spots on either side 
of the blue CO lobe (Fig. |b). This "high-velocity" CS(2- 
1) emission is likely to arise from dense, shocked material 
entrained by the outflow. 

The redshifted CS(2-1) emission detected on smaller 
scales by the PdBI interferometer in the 7-8 km ve- 
locity range also appears to be associated with the out- 
flow. This emission arises from the base of the red CO 
lobe and its shape closely follows the edges of the outflow 
lobe (see Fig. ^). No emission was detected by PdBI in 
the blue wing range (see Fig. ||d) or at the source sys- 
temic velocity. The blueshifted CS(2-1) emission seen in 
the single-dish map (Fig. ^) thus appears to be more ex- 



Based on the opacities derived above, optical depth effects 
cannot account for more than a factor of 1.3 decrease in the 
C^''S(2-1)/N2H+ (101-012) ratio. 
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3.6. Low degree of ionization 



Following, e.g., Caselli ct al. (1998), we can estimate 



the degree of ionization in the IRAM 04191 envelope 
from the observed value of the abundance ratio Rd = 
[DCO^]/[HCO''']. We measure an integrated intensity ra- 
tio /dco+/^hi3co+ = 3.3 ± 0.6 at the central position. 
Assuming optically thin R^^CO+{3-2) and DCO+(3-2) 
emission, this implies [DC0+]/[Hi3C0+] = 3.3 ± 0.6. 
Adopting a [^^C]/[^'^C] abundan ce ratio of 77 in the lo- 



cal ISM ( [Wilson fc Rood 1994| ), we thus derive Rd = 
0.04 ± 0.01, which is identical to the ratio measured by 
[Casclli et al. (200"lb ) toward the central position of L1544. 



As the same density (~ 10 cm~'^) is probed by the ob- 
servations in both cases, we conclude that the ionization 
degree Xi ^ 2 x 10~® derived in L1544 by Caselli ct al 



(2001b| ) based on their chemical models should be repre- 



sentative of the ionization degree in the IRAM 04191 en- 
velope at this density (assuming similar depletion factors 
for CO, which seems likely - see § |4.2| below). 



Fig. 6. CS(2-1) integrated intensity maps taken with the 4. Radiative transfer modeling: ID spherical 

30m telescope (left) and, on a smaller scale, with the PdBI simulations 

interferometer (right), (a) Single-dish map integrated in 

the 7.1 to 8.6 km s-^ velocity range (red), overlaid on Here, we use the radiative transfer code MAPYSO 



the C0(2-l) outflow contours from [AMB99| . (b) Same as 
(a) for the 4.6 to 6.1 km s~^ velocity range (blue). The 
CS wing emission (solid contours and gray scale) appears 
to be correlated with the CO outflow emission (dashed 
contours). The two crosses mark the central positions of 
IRAM 04191 [at (0,0)] and IRAS 04191. The CS(2-1) con- 
tour step is 0.3 K km s~^; the gray scale (shown on the 
right of each map) ranges from —0.3 to 3 K km s^^. (c) 
Cleaned interferometer map, integrated over the 7.1 to 8.6 
km s~^ velocity range (red), overlaid on the C0(2-l) con- 
tours of the outflow, (d) Same as (c) for the 4.9 to 6.1 
km velocity range (blue). The CS(2-1) contour step 
is 0.25 Jy beam~^ km s^^. The clean beam is shown on 
the bottom left. Note how the CS(2-1) red wing emission 
detected on a '--^ 5" scale coincides with the base of the 
red lobe of the CO outflow, while the CS(2-1) blue wing 
emission is almost entirely resolved out by the interferom- 
eter. 



tended than the redshifted emission and is likely resolved 
out by the interferometer. Such a difference in spatial ex- 
tent between blueshifted and redshifted emission, oppo- 
site to what infall motions would produce, is expected if 
the CS line wing emission arises from outflowing material. 
Likewise, the emission detected by the 30ni telescope near 
the systemic velocity, i.e., close to the dip of the CS(2-1) 
line (Fig. ^ , arises from extended foreground material on 
scales ~ 70" (see § ^3) and is also resolved out by PdBI. 



(Blinder 1997) to model the observed line spectra and set 
quantitative constraints on the kinematics of the IRAM 
04191 envelope (see Appendix for details about the code). 
For simplicity, we use a piecewise powerlaw description 
for the spatial variations of the kinetic temperature, den- 
sity, molecular abundance, infall velocity, and rotational 
velocity (in Sect. ||below) across the envelope (see Fig. 



4.1. Model inputs: mass distribution and kinetic 
temperature profile 



We use the env elope mass distrib ution derived by [Motte 
|fc Andre (200ll) (hereafter pAOl|) from 1.3 mm dust con- 
tinuum observations with the 30m telescope. MAOl esti- 
mate that the envelope mass contained within a radius 
of 4200 AU is Mo„v(r < 4200 AU) = 0.45 Mq with an 
uncertainty of a factor of 2 on either side, assuming a 
mean dust temperature Tjust — 12.5 ± 2.5 K obtained by 



AMB99 from a graybody flt to the A > 90 /im portion of 



the SED (see Fig. 3 of |AMB99| ). In addition, the radial 
structure analysis of MAOl indicates an average radial in- 
tensity profile / oc 0"™ with m ~ 0.6 ±0.1 in the range of 
angular radii 9 - 11" - 100" (i.e., r - 1500-14000 AU). 
Assuming a dust temperature profile Tdust oc r"* with 
q = -0.2 ± 0.2, |MA01| obtain a density profile p cx r-P 
with p = m+ l — 1.8 ± 0.3. Given the low bolomet- 
ric luminosity of IRAM 04191, the dust temperature is 
indeed expected to rise outward [q < 0) due to external 
heating by the interstellar radiation field (e.g. Masunaga] 
fe Inutsuka 200C; Evans ct al. 2001; Zucconi, Walmsley 



fe Galli 2001). The regime of central heating by the ac- 
creting protostar (q ~ 0.4) is likely confined to the inner 
~ 2000 AU radius region (cf. MAOl). Since a good thermal 
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Fig. 7. Input physical properties of the ID spherical en- 
velope model used to fit the observed CS and C'^'^S spec- 
tra (see Fig. P): Plots of density (a), column density (b), 
CS/H2 relative abundance (c), and kinetic temperature 
(d) as a function of distance from envelope center. 



coupling between gas and dust grains is expected for den- 
sities riHa > 10^ cm~'^ (e.g. Ceccarelli et al. 1996; Doty fc 



Neufeld 1997), the gas kinetic temperature profile should 
track the dust temperature profile, at least up to a radius 
of - 4000 AU. 

We also have some constraints on the gas kinetic tem- 
perature from our C^^O, N2H+, and CS observations. We 
measure a C^®0(2-1) to C^^O(l-O) integrated intensity ra- 



tio Ir 



C^"O(2-l)/^C^"O(l-0) 



= 1.5 ± 0.3 in the range of angular 
radii 0" — 60", using main beam temperatures and after 
degrading the resolution of the J = 2 — 1 data to that of 
the J = 1 — data. As the critical densities of C^*0(1- 
0) and 0^*0(2-1) are 3 x 10^ cm'^ and 2 x lO'' cm^^ 
at 10 K, respectively, both lines should be thermalized 
in most of the envelope. Assuming local thermodynamic 
equilibrium (LTE) and optically thin emission, we derive 
an excitation temperature Tex(C^*0) = 10 ± 2 K, which 
should be a good estimate of the gas kinetic tempera- 
ture in the low-density outer (r ~ 6000 AU) part of the 
envelope probed by C^^O. Likewise, the excitation tem- 
perature of the dense-gas tracer N2H+ (101-012) may be 
estimated from the relative intensities of the seven com- 
ponents of the N2H+(l-0) multiplet, assuming the same 
excitation temperature for all components (cf. § 3.3). This 
method yields Te,(N2H+) = 5.5 ± 0.5 K. With a criti- 
cal density '--^ 2 x 10^ cm"-^, the N2H+(l-0) multiplet is 
probably thermalized only in the inner (r < 2000 AU) 
envelope. We thus obtain a lower limit of ~ 6 K for the 
gas kinetic temperature in the dense (nn, > 10^ cm^"^), 
inner part of the envelope. Finally, the weak intensities 
of the optically thick CS(2-1) and CS(3-2) lines require a 



low gas kinetic temperature ~ 6-7 K in the range of radii 
2000-6000 AU. Given the density profile shown in Fig. ^i, 
a uniform gas temperature of 10 K would produce CS(2-1) 
and CS(3-2) spectra with main beam temperatures about 
1-2 K stronger than the observed temperatures. Beyond 
r ~ 6000 AU, the gas temperature is likely to increase to 
the typical ~ 10 K temperature of the Taurus cloud (see, 
e.g., Benson fc Myers 19891 ). However, the gas tempera- 
ture profile in the outer parts of the envelope has little 
influence on the CS and C'^'^S spectra since the observed 
lines are far from LTE there. 

Li summary, the gas kinetic temperature profile is 
likely to present a minimum of ~ 6 K at ~ 10" — 20" 
(i.e., r - 1400-2800 AU) and to reach a value of ~ 10 K 
in the outer parts of the envelope. 



4.2. Molecular depletion 

Assuming a standard isotopic ratio Xcs/Xc^^s = 22, a 
[CS]/[H2] abundance ratio of 8 x 10^^ is required to match 
the C'^'*S(2-1) integrated intensity at an angular radius 
40 " — 50 ". But a uniform abundance with such a value 
produces too strong G^^S{2-1), G^^S{i-2), and CS(5-4) 
spectra toward the center. A good fit to the C'^''S(2-1) 
integrated line intensities is obtained by assuming that the 
relative CS abundance drops by a factor of ~ 20 toward 



the center (see also § 3.4). Such a depletion factor for CS 
is comparable to those observed in starless cores such as 
L1544 (e.g. [Tafalla et al. 2002| ). 



4.3. Two regimes of in fa 1 1 

In Fig. H, we present a series of synthetic spectra emit- 
ted by a spherically symmetric model envelope with the 
input structure shown in Fig. ^ and our best estimate of 
the infall velocity field (shown by the solid line in Fig. [l^ a 
below). The model spectra are overlaid on the multitran- 
sition CS and C'^'^S spectra observed along the direction 
perpendicular to the outflow axis, so as to minimize the 
effects of the outflow. The blue infall asymmetry of the 
model optically thick lines and the position of the CS(2-1) 
dip match the observations well. The widths of the opti- 
cally thin lines are also well reproduced. The main short- 
coming of the model is that it does not reproduce the fairly 
strong emission present in the wings of the observed CS 
spectra. As discussed in § |3.5| , we attribute these wings to 
the fraction of envelope material entrained by the outflow. 

In order to determine the range of input model pa- 
rameters that yield reasonably good fits to the observed 
CS and C'^^S spectra, we have performed a comprehen- 
sive exploration of the parameter space, as illustrated in 
Figs. ||, |l^ and |ll|. Strong constraints on the infall velocity 
arise from the small optical depth lines, i.e., C'^'*S(2-1) and 
C34S(3-2). The widths of these hues set firm upper limits 
to the absolute value of the infall velocity on the size scale 
of the beam. We obtain Vini < 0.15 km s~^ at ~ 1750 
AU and Vini < 0.2 km s"^ at ~ 1150 AU from C34S(2-1) 
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Fig. 8. CS(2-1), CS(3-2), CS(5-4), C34S(2-1), and C34S(3-2) spectra (in units of main beam temperature) observed 
with the 30m telescope along the direction perpendicular to the outflow axis (histograms). The dotted line indicates 
our best-fit estimate (6.63 km s^^) of the envelope systemic velocity based on o ur CS/C'^'^S modeling. Synthetic 



spectra corresponding to the 'best-fit' 
model parameters) are superposed. 



ID spherical collapse model described in § 13 (cf. Fig. and Fig. [l^a,b for 
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Fig. 9. Influence of the outer infall velocity field on the CS(2-1) {upper row) and C'^'*S(2-1) {lower row) line profiles. 
Synthetic spectra (solid curves) corresponding to a range of model infall velocities are superposed on the spectra 
observed at the central position of IRAM 04191 (histograms). All models have Vmf = vq = 0.2 km s~^ at ro = 750 AU, 
Vinf oc r~°'^ for r < rbreak, and Vinf = Vcruise for r > rtreak- The models differ by the value of rbreak, or equivalently by 



the value of the infall velocity at large radii, Vcmise = v^y. {—^j '■ rbreak ~ 12000 AU (i.e., Vcruise ~ 0.05 km s ) 

in (a), rbreak — 3000 AU (i.e., Vcruise = 0.1 km s^^) in (b) (= preferred model shown in Fig. ||), and rbreak = 750 
AU (i.e., Vcruise = 0.2 km s~^) in (c). The dotted vertical line marks the envelope systemic velocity as in Fig. ||, 
while the dashed and dash-dotted lines mark the velocities of the dips in the observed and synthetic CS(2-1) spectra, 
respectively. Note that the outer infall velocity Vcmise is too small in model (a) and too large in model (c) to match 
the observed redshift of the CS(2-1) dip and the width of the €^''8(2-1) hue. 



and C^"'S(3-2), respectively. On the other hand, the am- 
plitude of the blue asymmetry seen in the self-absorbed 
CS(2-1) and CS(3-2) lines, as well as the redshifted posi- 
tion of the corresponding absorption dips (see § 3.3), both 
require a relatively flat, extended infall velocity field with 
^;inf ~ 0.10 ± 0.05 km s'^ up to r - 10000 - 12000 AU. 
The latter value approximately corresponds to the radius 
where the bulk of the absorption occurs in CS(2-1) and 
CS(3-2) (see § p^ ). These constraints are illustrated in 
Fig. H which shows the effect of varying the infall velocity 
field on the central CS(2-1) and 0^48(2-1) spectra. Three 
models are compared: the preferred model displayed in 
Fig. H is shown in the central panel (Fig. |9|3) , while models 
with lower and higher infall velocities in the outer part of 
the envelope are shown in the left (a) and right (c) pan- 
els, respectively. It can be seen that the position of the 
CS(2-1) absorption dip is not redshifted enough in model 
(a) and too redshifted in model (c) to match the obser- 
vations. Furthermore, the CS(2-1) blue-to-red asymmetry 
is too weak in model (a) and the C'^'*S(2-1) line becomes 
too broad in model (c) compared to the observations. Only 



model (b) approximately reproduces the observed position 
of the CS(2-1) dip and the width of the C34S(2-1) line. 

In the context of a pure infall model, larger veloci- 
ties (winf ~ 0.2 - 0.4 km s"^) in the inner {r ~ 700 AU) 
part of the envelope are suggested by the broad linewidth 
of the central CS(5-4) spectrum. This is shown in Fig. [l^ 
which compares three models differing in the magnitude of 
their infall velocity at 700 AU. A similar trend is indicated 
by the broadening of the N2H+ (101-012) line toward the 



center (cf. § 3.4). However, the central CS(5-4) spectrum 
(which does show infall asymmetry - cf. Fig. ^) may be 
partly contaminated by small-scale structure in the out- 
flow. In particular, the redshifted portion of the CS(5-4) 
spectrum, not reproduced by the model of Fig. ^, may be 
related to the redshifted CS(2-1) emission detected by the 
interferometer on small scales (see § 3.5). Such a switch 
between infall-dominated CS emission on large scales and 
outflow-dominated CS emission on small scales is also ob- 



served in the Class object B335 (cf. |Wilncr et al. 2000|) 



In addition to infall, a "turbulent" velocity field is 
needed to match the width of the dip, and consequently 
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Fig. 10. Influence of the inner infall velocity field on the CS(5-4) {upper row) and C'^'*S(2-1) {lower row) line profiles. 
Synthetic spectra (solid curves) corresponding to a range of model infall velocities are superposed on the spectra 
observed at the central position of IRAM 04191 (histograms). All models have Vinf esc r~°-^ for r < r^reak and Vinf = 0.1 
km s~^ for r > rbreak- The models differ by the value of rbreak, and therefore by the value of the infall velocity at 
the 700 AU radius of the CS(5-4) beam: rtreak = AU (i.e., wtooau = 0.1 km s~^) in (a), rhreak — 3000 AU (i.e., 
?^700AU = 0.2 km s~^) in (b) (= preferred model shown in Fig. ||), and rhreak = 13900 AU (i.e., W700AU — 0.4 km s~^) 
in (c). The dotted vertical line marks the envelope systemic velocity as in Fig. |[ Note that the CS(5-4) linewidth is 
too small in model (a) and the C'^^S(2-1) linewidth is too large in model (c), suggesting 0.1 < Vinf < 0.4 km s~^ at 
r - 700 AU. 
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Fig. 11. Influence of the turbulent velocity field on the CS(2-1) spectra. The line profiles observed at the central 
position of IRAM 04191 (histograms) are compared to synthetic spectra (solid lines) corresponding to models with a 
uniform turbulent velocity for various values of the velocity dispersion oturh'- (a) Oturh = km s~^, (b) aturb = 0.085 
km s^^ (= reference model shown in Fig. |^), (c) aturb = 0.13 km s^^. The dotted vertical line marks the envelope 
systemic velocity as in Fig. |[ Note that the dip of the model CS(2-1) spectrum is too narrow in (a) and too broad in 
(c). 
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the velocity difference between the blue and the red 
peak, in the optically thick CS(2-1) and CS(3-2) spec- 
tra. Assuming a uniform turbulent velocity dispersion for 
simplicity, a good compromise between the upper limit 
set by the linewidth of the optically thin C'^^S(2-1) and 
C^'*S(3-2) spectra and the lower limit set by the width 
of the CS(2-1) and CS(3-2) dips is obtained for crturb = 
0.085 ± 0.02 km s^^ (cf. F ig. |ll|). This is equivalent to 
AwfW,HM = (Tturb X V81n2 = 0.20 ± 0.05 km s'^ and cor- 
responds to only half the thermal broadening of the mean 
molecular particle at 10 K, showing that the IRAM 04191 
envelope is "thermally-dominated" (see also § 3.4) as are 



Taurus dense cores in general (e.g. Myers 1999). 

The main conclusions of our ID exploration of the pa- 
rameter space are summarized in Fig. [12^ and Fig. p^b, 
where the shaded areas represent the ranges of infall ve- 
locities (a) and turbulent velocity dispersion (b) for which 
acceptable fits are found. Two infall regimes seem to 
stand out in Fig. Oa: the infall velocity is relatively large 
('^inf ~ 0.2 km s~, supersonic) and increases toward the 
center for r ;$ 2000 — 3000 AU, while it is smaller and 
roughly uniform at Winf '^-^ 0.10 ± 0.05 km s~^ between 
- 2000 - 3000 AU and - 10000 - 12000 AU. Given the 
density profile of Fig. such an infall velocity field im- 
plies a mass infall rate of Mini '--^ 3 x 10~^ M© yr~^ at 
r = 1750 AU. (The density and velocity profiles shown in 
Fig. ^ and Fig. ^2|a are such that Afinf is roughly inde- 
pendent of radius.) Inside the r ~ 11000 AU region (where 
non-zero inward motions are inferred), the fraction of en- 
velope mass with supersonic ( ^ 0.16 — 0.2 km s^^) infall 
velocities is estimated to be only ~ 1 — 10%, depending 
on the exact value of the sound speed and exact form of 
the infall profile (see Fig. [T^a). 



5. Radiative transfer modeling: Simulations with 
infall and rotation 

5.1. Quasi 2D simulations 

T o ac count for the effects of rotation in the envelope (see 
§ and Fig. ||), we have performed "quasi" -2D simu- 
lations with the following approximation. The non-LTE 
level populations arc still calculated with a ID Monte 
Carlo method (see Appendix) assuming a spherical enve- 
lope with the same characteristics as the model described 
in § 4.3. We then add a cy lindrical rotation velocity field 
to the ID model of S 4.3 and use the 2D version of the 



MAPYSO code to compute a proper radiative transfer in- 
tegration along each line of sight. If we ignore departures 
from a spherical density distribution, this approach would 
remain strictly exact in the case of solid-body rotation, 
since the velocity difference between any couple of points 
projected on the axis joining these points is insensitive 
to the addition of a solid-body rotation component (cf. 
Ward-Thompson & Buckley 2001). In practice, however, 
the rotation observed here departs from solid body and 
the density distribution is not spherical. We therefore as- 
sume that, to first order, the line excitation is much more 
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Fig. 13. (a) C'^''S(2-1) position-velocity diagram of the 
'best-fit' 2D collapse model including differential rota- 
tion, taken along the direction perpendicular to the ro- 
tation axis, (b) Observed C'^^S(2-1) position-velocity di- 
agram along the direction perpendicular to the outflow 
axis. In both diagrams, the dots represent the peak veloc- 
ities derived from Gaussian flts at angular radii larger than 
20 " and from centroid (flrst moment) estimates for radii 
smaller than 20 "; the curve shows the variation of the 
peak velocity (estimated in the same way) in the model. 

sensitive to the density distribution (averaged over angles) 
than to the velocity field. (In particular, we have checked 
that the profiles of excitation temperature are essentially 
insensitive to the infall velocity field as long as the turbu- 
lent velocity dispersion is of the same order as the infall 
velocity.) The rotation velocity field has nevertheless im- 
portant effects on the shape of the line profiles, which we 
properly take into account here. 



5.2. Two regimes of rotation 

The cylindrical rotation velocity field that we have added 
to the ID spherical model of 



4.3 is shown in Fig. \L 



The rotation axis is taken to coincide with the outflow 
axis, at a n inclination angle i = 50° to the line of sight 
(see § L2) and a position angle P. A. = 28 ° in projection 
onto the plane of the sky (cf. Fig. |l|). 

Adding rotation yields two major improvements in the 
fits to the CS and C^^S spectra (see Fig. |l| and Fig. |l|). 
First, the velocity gradient of the small optical depth 
C'^^S(2-1) line along the direction perpendicular to the 
outflow is well reproduced, as shown in Fig. ^ Both the 
model C^"'S(2-1) position-velocity diagram (Fig. 13a,) and 
the centroid velocity curve agree well with the observa- 
tions (Fig. 13 b). 

Second, the blue asymmetry of the model CS(2-1) and 
CS(3-2) spectra is enhanced toward the north-west and 
attenuated to ward the south-east, as seen on the observed 
spectra (see § 3^ and Fig. |lj). 

In agreement with the discussion of the position- 



velocity diagrams (§ 3.2 above), the present 2D modeling 
indicates that the envelope can be divided into two regions 
with distinct rotational characteristics (see the radial pro- 
files derived for the rotational velocity and angular veloc- 
ity in Fig. |2fc and Fig. p^, respectively). First, solid-body 
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Fig. 12. Infall (a), turbulence (b), and rotation (c) velocity fields inferred in the IRAM 04191 envelope based on our 
ID (§ ^) and 2D (§ 0) radiative transfer modeling. The shaded areas show the estimated domains where the models 
match the CS and Cr^'S observations reasonably well. In (a) and (b), the solid lines show the infall velocity and the 
turbulent velocity dispersion in both the ID and 2D models (cf. Fig. || and Fig. |l^, respectively) as a function of radius 
from envelope center. In (c), the solid line represents the profile of the azimuthal rotation velocity in the 2D envelope 
model (cf. Fig. |l|) as a function of distance from the outflow/rotation axis. The point with error bar at 11000 AU 
corresponds to the velocity gradient observed in C^^O (cf. § |3.2| ). Panel (d) shows the corresponding angular velocity 
profile. 



Vrot oc r 

-2.5±0.5 



rotation is ruled out in the outer 3500 < r < 7000 AU ra- 
dius envelope, where a good fit to the centroid velocity 
curve is obtained with v- 
an angular velocity f2 oc r 

r < 3500 AU radius region, our simulations suggest a ve- 
locity profile Urot 7-0 i±0-4^ j fi qt ^-o.9±o.4^ although 
the form of the position-velocity diagram is significantly 
influenced by the finite resolution of the observations. A 
rotation velocity Wi-ot — 0.20 ±0.04 km s^^, corresponding 



i.5±o.5^ corresponding to 
Second, in the inner 
ions sug 

-0.9±0.4 



to = 12±3 km s pc , is derived at a radius r = 3500 
AU after correcting for inclination^ 



Tiie reason why the value of 12 km s ^ pc ^ found here for 
Q, is larger than the apparent velocity gradient of 9 km pc" 



given in § 3.2 is that the 2D model shown in Fig. |l^ properly 
accounts for the finite (~ 25") resolution of the observations. 



Belloche et al.: Velocity structure of the IRAM 04191 envelope 



15 




Fig. 14. Same as Fig. p| b ut with synthetic spectra corresponding to the 'best-fit' 2D model including differential 
rotation described in § ^.21 (cf. Fig. ^ and Fig. ^ for model parameters). 
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6. Discussion: Comparison with collapse models 



ence (e.g. Tomisaka 1996; Basu 1997; safier, McKee & 
Stabler 1997| ; |Li 1998| ; |Ciolek fc Konigl 1998D of magnetic 



In this section, we fli'sL sunmiai'ize Llie pi'cdlcLiouy of ml- 
lapse models and the main constraints derived from our 

|) , and then disc uss t he appli- 

~|to§ii). 



observations 



6.1 and 

cability of various models to IRAM 04191 
Some implications for the distribution and evolution of 
angular momentum during protostellar collapse are dis- 
cussed in S Bii. 



6.1. Overview of model predictions 

Whitworth & Summers (1985) have shown that there is a 
two-dimensional continuum of similarity solutions to the 
problem of isothermal spherical collapse. In this contin- 
uum, the well-known solutions proposed by ^hu (19"77 ) 
and Larson (1969 )- Pension (1969 ) represent two extreme 
limits. All isothermal similarity solutions share a univer- 
sal evolutionary pattern. At early times (i < 0), a com- 
pression wave (initiated by, e.g., an external disturbance) 
propagates inward at the sound speed, Cs, leaving behind 
it a p{r) oc density profile and a uniform infall ve- 
locity field. This compression wave has zero amplitude in 
the limiting case of the Shu 'inside-out' collapse solution. 
At t = 0, the compression wave reaches the center and 
a point mass forms which subsequently grows by accre- 
tion. At later times {t > 0), this wave is reflected into a 



rarefaqtion or expansion wave, propagating outward (also 
at the sound speed) through the mfallmg gas, and leav- 
ing behind it free-fall density and velocity distributions 
(i.e., p(r) oc r~^'^ and v{r) cx r~^'^). The various solu- 
tions can be distinguished by th e absolute values of the 
density and velocity at t 0. The Shu (1977 ) solution has 
the p{r) = ifls'^ /2tiG) density distribution of a static 
{v = 0) singular isothermal sphere (SIS) at t = 0, while 
the Larson-Penston (1969) solution is ~ 4.4 times denser 
and far from equilibrium (w « —3.3 Os). The recent find- 
ing of inward motions of subsonic amplitude ~ 0.02 — 0.10 
km s~\ extended over ~ 0.1 pc (~ 20000 AU) in the 



prestellar core L1544 (see Tafalla et al. 1998 ; Williams fc 
Myers 1999 ; [Caselli et al. 200"Ta ) suggests that true pro 



tostellar collapse in Taurus proceeds in a manner which is 
neither the Shu nor the Larson-Penston flow, and is per- 
haps more reminiscent of an intermediate similarity solu- 
tion. 

In practice, the initial conditions for fast protostellar 
collapse are not strictly self-similar and involve a density 
profile that is flat at small radii (e.g. Ward-Thompson 
et al. 1994| ; [Andre, Ward-Thompson fc Motte 1996] ) and 
bounded or sharp-edged at some outer radius i?out (e. 



Motte, Andre & Neri 1998 




Bacmann et al. 2000 


) like a 


finite-sized Bonnor-Ebert isothermal sphere (e.g. 


Bonnor 



195£; Alves et al. 2001). A number of recent numeri- 
cal (magneto)hydrodynamic simulations or simplified an- 
alytical calculations attempt to describe the collapse in 



such a situation, either in the absence (e.g. 


Foster & 


Chevaher 1993 


; Henriksen et al. 1997 


; Masunaga, Miyama 


& Inutsuka 199? ; 


Hennebelle et al. 2002 


) or in the pres- 



fields. The Larson-Penston similarity solution is found to 
describe the collapse quite satisfactorily near i = (at 
least for small radii), but the Shu solution is more ade- 
quate at intermediate t > times, before the expansion 
wave reaches the edge of the initial, pre-coUapse dense 
core. 

When rotation is included, a rotationally-supported 
disk develops at the center of the infalling envelope during 
the accretion phase (i.e. at t > 0). The size scale of this 
disk is determined by the centrifugal radius, Rc, which 
defines the position where the centrifugal force balances 
gravity in the equatorial plane. Strong departures from a 
spherical density distribution and a purely radial inflow in 
the envelope are expected to occur on size scales of order 
(or smaller than) Rc (cf. [Chevalier 1983 ; Hartmann 1998). 
Most collapse models predict that Rc = jGm should in- 
crease with time as material of higher and higher specific 
angular momentum j = fl R^ falls in, but the exact de- 
pendence on time i, or alternatively accumulated central 
mass TO, varies from model to model, according to the dis- 
tributions of mass and angular momentum at t = 0. For 



instance, Rc scales as or in the Shu model (TSC84), 
which assumes solid-body rotation at point mass forma- 
tion. By contrast, the dependence of Rc on to is only lin- 
ear in the magn etically-controlled collapse mo del of Basi; 
(1998[ ) (see also [Krasnopolsky fc Konigl 2002[) . 



6.2. Summary of observational constraints 

The analysis of our line observations (§ ^ and § ||) indi- 
cates that both the infall and rotation velocity fields of 
the IRAM 04191 envelope are characterized by an inner 
and an outer regime (see Fig. |2|). 

Our ID radiative transfer simulations (§ ^ indicate that 
the infall velocity profile is flat with Vi^f ~ 0.1 km s^^ be- 
tween n - 2000 AU and n^o 10000 - 12000 AU. Higher 
infall velocities at radii r < ri are suggested the CS obser- 
vations, which are consistent with a free-fall velocity field 
{vinf OC r~°-^) at r < Ti. The width of the optically thin 
C'^^'S lines strongly constrains Vini to be ^ 0.15 km s^^ 
at r - n (cf. Fig.|l2|a). 

The position-velocity diagrams observed in optically 
thin lines (§ 3.2) show that the envelope is differentially 
rotating with an angular velocity profile (x 7--2.5±o.5 
tween r,n 3500 AU and r^.o 7000 AU, and ~ 12 
km s~^ pc~^ at r ~ r^^. Although the limited spatial 
resolution of our observations prevents us from deriving 
accurate values in the inner region, the rotation profile is 
definitely shallower for r < rm (cf. Fig. [l^). A follow-up 
interferometric study is underway to provide more accu- 
rate constraints on the angular velocity in this inner re- 
gion. 

Although the two critical radii and r„i differ by less 
than a factor of ~ 2, it is unclear whether they are physi- 
cally related or not. The inner r < rj envelope may corre- 
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Fig. 15. Ratio of centrifugal to gravitational acceleration 
inferred in the IRAM 04191 envelope based on our 2D 
(§ H) radiative transfer modeling. The dashed area shows 
the estimated domain where the models match the CS and 
C^'^S observations reasonably well. The solid line shows 
the pr ofile of amntr I g,- in the 2D envelope model presented 



in Fig. 12, 13 and 14 , assuming ^ 0.05 Mg,. 



radius = 3500 AU of the rapidly rotating inner enve- 
lope, showing that the centrifugal acceleration is a sizeable 
fraction of the gravitational acceleration in this region. 
Comparable values of ftcontr/ffr are nevertheless obtained 
in some magnetic models of cloud collapse (see Fig. 4 of 
Basu fc Mouschovias 199^ and § |6.5| below) . 

We estimate the centrifugal radius to be Rc < 400 AU, 
assuming the maximum rotation velocity profile consis- 
tent with the observations (cf. Fig. [l^c) and a stellar 
mass Mi, = 0.03 M©. In principle, a large centrifugal ra- 
dius should imply the pres ence of a large centri fugal disk 
around the central object, ^tahler et al. (1994 ) have in- 
vestigated the initial growth of protostellar disks in the 
context of the rotating collapse picture of TSC84 and have 
shown that a 3-component structure should develop inside 
Rc'. a dense inner accretion disk of radius i?disk ~ 0.3iRc, 
a ring where material and angular momentum pile up at 
i?disk, and a low-density outer disk where material travels 
at high velocity between Rc and i?disk- If -Rc ~ 200 AU, we 
may thus expect i?disk ~ 70 AU. Our 1.3 mm continuum 
interferometric observations set a firm upper limit to the 
radius of any dense inner circumstellar structure around 
IRAM 04191: iJdisk < 10 AU (cf. § |]l]). This suggests 
that the actual centrifugal radius is significantly smaller 
than 200 AU or else that the presence of a protostellar 
companion at r ~ 30 AU may have cleared a substan- 
tial gap in th e inner accretion disk (cf. [lensen, Mathieu fc 
Fuller 19961). Artymowicz & Lubow (1994) show that, due 



spond to the free-fall region developing inside the expan- 
sion wavefront, while the outer r > rj envelope may be 
dominated by the flat, extended inward velocity field set 
up by the compression wave at < < (see § 3.1 above). If 
the current radius of the expansion wave is indeed ~ ri, 
then the age of the central protostellar object should be 
t ^ n/fls 5 X 10"* yr (assuming a propagation speed 
Us ~ 0.2 km s^^), in r ough agreement with the estimated 
age of ;S 3 X 10** yr (§ [L^ ). On the other hand, the fastly 
rotating r < r,^ region may be a dynamically collapsing 
'supercritical' core in the process of decoupling from the 
ambient medium, and the outer < r < rm.o envelope 
may be a transition regi on b etween the protostar and the 
background cloud (see § 6.5 below). 



to tidal disk truncation in binary systems, any individual 
disk must have an outer radius less than half the binary 
separation, and any circumbinary disk must have an inner 
radius more than roughly twice the binary separation. The 
non-detection of a dense ~ 70 AU radius disk by the PdBI 
interferometer could thus be accounted for by this effect 
if IRAM 04191 were a protobinary of separation ~ 30 AU 
(i.e. - 0.2"). 



6.3. Problems with the inside-out collapse model 

The inside- out, isothermal collapse model described by 
^hu (19*77 ), and its 2D extension including rotation 



(TSC84), has been widely used to explain the infall spec- 



tral signature seen in the envelopes of low-mass protostars 



(e.g. Zhou et al. 1993; 


Myers et al. 199£; Choi et al. 1995 


^houetal. 1996; 


Hogcrheijde & Sandell 200C). It accounts 



6.2.1. Centrifugal support 

Using our derived model for the structure and kinematics 
of the IRAM 04191 envelope (e.g. Fig. and Fig. |l|), 
we can estimate the dynamical importance of rotation 
in the envelope. The ratio of the centrifugal accelera- 
tion acontr = "^lotl"^ to the local gravitational field = 
G X [Menv('") M^]/r2 is shown in Fig. |l|. Here, Me„v(r) 
is the envelope mass within radius r corresponding to the 
density profile of Fig. and Af^ is the uncertain mass 
of the central protostellar object ^ 0.03 — 0.1 Mq (cf. 



relatively well for the densities and accretion rates mea- 



sured in Taurus protostellar envelopes (e.g. Ohashi 1999 



Motte & Andre 2001; this paper). This model is, how- 
ever, inconsistent with the combined density and infall 
velocity profiles measured for IRAM 04191 (cf. Fig. |^a 
and Fig. ^2|a). An expansion wave radius <L 10000 AU 
is indeed required to reproduce the blue asymmetry of 
the central CS(2-1) spectrum in the context of the inside- 
out collapse picture. But such a large infall radius then 
implies high infall velocities at small radii r ^ 1500 AU 
which in turn yield central 03*8(2-1) and 03*8(3-2) spec- 
AMB99). We estimate a ratio flcontr/Sr ~ 0.4 ± 0.2 at the tra that are 2-3 times too wide compared to the observed 
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linewidths. Conversely, a model with an expansion wave 
radius of only ~ 2000 AU yields correct C34S(2-1) and 
C^'*S(3-2) linewidths but fails to reproduce the strong 
asymmetry observed in CS(2-1) and CS(3-2). It is clearly 
because of the absence of significant inward motions at 
t < that the inside-out collapse model cannot account for 
the infall velocity profile observed here. A similar incon- 
sistency with the inside-out collapse model has been noted 



much better fit to the observations than either the Shu or 
the Larson-Penston similarity solution. A definite problem 
with the critical Bonnor-Ebert model, however, is that as 
much as 44% of the envelope mass is predicted to flow in 
supersonically at t = (Foster & Chevalier 1993), and pro- 
gressively more mass at t > 0. This is much larger than the 
1 - 10% mass fraction derived in S Ofor IRAM 04191. 



previously in prestellar cores such as LI 544 by Tafalla et 
al. (19f|) and |Lee, Myers fc Tafalla (200l| ). 



6.4. Comparison with other thermal models of collapse 

With an age < 3 - 5 x 10"* yr for the central IRAM 04191 
protostellar object (see § 1.2 and § 6.2) and a sound speed 



Qs ~ 0.19 km s~^ (Tk — 10 K), the radius of the expansion 
wave must be smaller than the inner radius rj ~ 2000 AU 
of the extended region where a flat infall velocity field is 
inferred (see Fig. [l2|a). The kinematics and density struc- 



ture o | the outer r > n envelope should thus still re- 
flect the physical conditions at t = 0. Accordingly, the 
observation of substantial infall velocities between ri and 
ri_o ~ 11000 AU points to collapse models that are more 
dynamical than the Shu solution and involve the propa- 
gation of a finite-amplitude compression wave prior to the 
formation of the central object (cf. § S.l and Whitworth 



We conclude that spherically symmetric collapse models 
in which thermal pressure provides the only force opposing 
gravity tend to be too dynamical and are only marginally 
consistent with our observational constraints. It is likely 
that the inclusion of rotation in thermal models would 
improve the comparison with observations and may even 
account for the flattened shape of the protostellar envelope 
seen in Fig. [| (P. Hennebelle, private communication). 
However, a generic feature of rotating thermal collapse 
models is that, due to simultaneous conservation of en- 
ergy and angular momentum, they tend to predict similar 
forms for the rotation and infall velocity profiles at radii 
(much) larger than the centrifugal radius (e.g. Baigo fc 
Hanawa 199^ ) . This is at variance with the steeply declin- 



& Summers 1985). On the other hand, the infall velocities 
derived in the outer (r > n) envelope are only subsonic 
(approximately half the sound speed) , and clearly incon- 
sistent with the supersonic infall velocities ~ 3.3aa charac- 
terizing the Larson-Penston isothermal similarity solution 
at i = (cf. § 6.1). The infall velocity field derived here 
(cf. Fig. [T^a) is suggestive of a moderately dynamical col- 
lapse model, intermediate between the Shu similarity so- 
lution (zero-amplitude compression wave at i < 0) and 
the Larson-Penston solution (strong compression wave at 
t < 0). 

Qualitatively at least, such a moderately dynamical in- 
fall velocity field resembles that achieved during the col- 
lapse of a finite-sized, Bonnor-Ebert isothermal sphere. 
For instance, in their numerical simulations of the collapse 
of critically stable Bon nor-Ebert spheres (withou t mag- 
netic field or rotation), Foster fc Chevalier (1993 ) found 
infall velocities at i = ranging from 3.3as near the ori- 
gin to at the outer boundary radius. Assuming that the 
IRAM 04191 dense core was initially a marginally stable 
Bonnor-Ebert sphere with a center-to-edge density con- 
trast of ~ 14, the central density at the onset of collapse 
must have been ~ 5 x 10"* cm"'^, given the typical outer 
density ^ 3 x 10"^ cm^^^ observed in prestellar cores (e.g. 



ing rotation velocity profile and flat infall velocity profile 
we observe beyond 3500 AU (a radius much larger than 
the estimated centrifugal radius - see § |6.2.1 above). In 
fact, the strong decline of the rotation velocity profile be- 
yond 3500 A U suggests that angular momentum is not 
conserved in the outer envelope. In the next subsection, 
we propose that this results from magnetic braking. 



6.5. Comparison with magnetically-controlled collapse 
models 



Magnetic ambipol a r diffusion models 
Mouschovias 1994: Basu 



Ciolek 



Mouschovias 1994, 1995a 



- hereafter [BM94 |BM954 |BM95bD are another class of 
models which yield non-zero inward velocities in an ex- 
tended region prior to point mass formation. Am bipolar 
diffusion has been invoked by Ciolek fc Basu (200C ) to ex- 
plain the extended inward motions observed in the Taurus 
prestellar core L1544 (see Tafalla et al. 1998 ; Williams fc 
[Myers 1999| ; [Caselh et al. 2001aD . The models start from a 



Bacmann et al. 2000). Adopting a temperature of 10 K, 
Fig. 1 of Foster fc Chevalier (1993) then predicts an infall 



velocity aX t — Q varying from ^ 0.2 km s at ~ 6000 



AU to 0.5o 



0.1 km s- 



at 



15000 AU. Even if this 



represents a somewhat stronger variation of infall veloc- 
ity with radius than derived in the case of IRAM 04191 



(Fig. [13a), the Foster & Chevalier model clearly provides a 



magnetically subcritical cloud, initially supported against 
gravitational collapse by a static magnetic field, and pre- 
dict an evolution in two phases. During the first, qua- 
sistatic phase, the subcritical cloud contracts along direc- 
tions perpendicular to the field lines through ambipolar 
diffusion. The gravitationally-induced inward drift of the 
neutral species is slowed down by collisions with the ions 
which are well coupled to the magnetic field. The cen- 
tral mass to magnetic flux ratio increases with time, un- 
til it reaches the critical value for collapse, {l/2^T)G~^^^ . 
The inner region of the cloud then becomes magnetically 
supercritical and collapses dynamically, while the outer 
envelope remains subcritical. This effect introduces a spa- 
tial scale in the collapse process, which corresponds to 
the boundary between the magnetically supercritical in- 
ner core and the subcritical outer envelope. These two re- 
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gions are characterized by distinct rotational properties. 
The supercritical inner core evolves with conservation of 
angular momentum and rapidly spins up. It achieves a 
power-law angular velocity profile ^(r) oc S(r) oc 1/r at 
t = 0, where E is the (mass) column density (e.g. Basu 



1997|)~|By contrast, due to magnetic braking, the subcriti- 
cal envelope loses (part of) its angular momentum on the 
timescale tj « iyr^ (where t^y is the free-fall time at the 
density of the background - cf. Spitzer 1978 and Tomisaka 
2000jy^nd is progressively brought to near corotation with 
the background medium, assumed to rotate at the uniform 
rate fib (e.g. BM94). This separation generates a break in 
the angular velocity profile at the radius Rcrit of the mag- 
netically supercritical core. 

Qualitatively at least, the rotation and infall profiles 
of the IRAM 04191 envelope (see Fig. Il| and § K% can 



be accounted for in the framework of such magnetic mod- 
els[[| if we identify Rcrit with the radius ~ 3500 AU be- 
yond which the observed rotation profile exhibits a marked 
steepening (cf. Fig. p^). Indeed, the angular velocity pro- 
file of the inner r < region is consistent with the 
17 oc 1/r power law expected in the supercritical core at 
t ~ 0. Furthermore, the models predict a steepening of the 
rotation profile beyond Rcrit when magnetic braking does 
not have enough time to bring the system to corotation 
with the background before the formation of the super- 
critical core. In the parameter study presented by BM95a 
and BM951: , this happens in models 6 and 8, for two differ- 
ent reasons. First, if the cloud is already critical (or close 
to critical) initially near the center (but still subcritical 
in its outer parts) as in model 6, there is no quasistatic 
phase and the dynamical contraction of the supercritical 
core starts right away on a timescale ^ iyr^^ (free-fall time 
at cloud center) shorter than tj. This produces a tran- 
sition region beyond Rcrit with a steep angular velocity 
profile between the supercritical inner core and the back- 



ground outer region (see Fig. 6b of BM95b). At the same 
time, the collapse of the supercritical core is retarded by 
the magnetic forces so that supersonic infall velocities de- 
velop only very close to the center (at r ^ 0.01 x Rcrit)- 
Second, if the ionization fraction is low enough that the 
magnetic braking timescale is only slightly shorter than 
the ambipolar diffusion timescale as in model 8 of BM95a, 



then the supercritical core can begin its dynamical evolu- 
tion with a rotational angular velocity much larger than 
that of the background. This also results in the appari- 
tion of a transition region with a steep rotation profile be- 
tween the sup ercritica l core and the e xternal b ackground 
(cf. Fig. 5 of |BM95a|). In model 8 of |BM95at the outer 



radius of the transition region is ^ 3 x Rcrit , which is con- 



Strictly speaking, the ambipolar diffusion models of Basu 
& Mouschovias describe the evolution of core properties only 
during the prestellar phase (t < 0), while the properties of the 
IRAM 04191 envelope beyond ri ~ 2000 A U ar e representative 
of the physical conditions at t ~ (cf § 6.4). However, the 
models are easily extrapolated to f = by considering th e lim- 
iting profiles achieved for an infinite central density (cf. Basrj 
199|, |199^). 



sistcnt with the extent of the oc r~^-^ zone in Fig. [ijd. 
In this model, the infall velocity becomes supersonic for 
r ^ 0.3 X Rcrit and exhibits a relatively fiat profile beyond 
Rcrit I E^lso in agreement with the observational constraints 
of Fig. ^2|a. The ionization degree Xi = ni/n(R2) assumed 
in the model is low (2 x 10~^ at a density of 10^ cm~'^), 
but comparable to the value derived at the center of the 
IRAM 04191 envelope in 



3.6 



Quantitatively, however, it is more difficult to obtain 
a good match of the observations of IRAM 04191 with 
published ambipolar diffusion models. These models ro- 
tate a factor of ~ 3 — 10 more slowly^ than does the 
IRAM 04191 envelope and have supercritical core radii 
that are a factor ~ 3 — 30 bigger than the observed 
break radius r,n ^ 3500 AU (cf. Table 2 of |BM95bD . 
Physically, the radius of the magnetically supercritical 
core corresponds to the Jeans length for the density Ucrit 
(or equivalently surface density T,crit) at which "decou- 
pling" occurs, i.e., the density (or surface density) at 
which gravity overcomes magnetic support. In the disk- 
like geometry of ambipolar diffusion core models, the crit- 
ical Jeans radius is Rcrit = 



ncrit — 2timH 

BM95a|) . For Tcore 
gives Ncrit = — 



X (Ecrit/as) 



= 7 K and Rcrit 
■ ~ 9.3 X 10^1 c 



I {2GY,crit) and one has 
X {aJRcritY (see 
^ 3500 AU, this 
m^^ and Ucrit ~ 



1.4 X 10^ cm~^, the latter being remarkably similar to the 
volume density estimated at r = 3500 AU in the envelop e 
(as expected in the models - see, e.g., Fig. 6a of BM95b ). 
If the mass-to-flux ratio Af/$ = E/B is just critical in 
the supercritical core, then the magnetic field strength 
of the core should be Bcrit = 27rG^/^Ecrit ~ 60 /iG at 
r 3500 AU. The reason why published magnetic mod- 
els have larger values of Rcrit is that their critical "de- 
coupling" densities and field strengths are typically lower 
than these estimates by factors ~ 10 and ^ 3, respec- 
tively. Direct Zeeman measurements suggest that the mag- 
netic field strength is only ~ 10 /j,G in the low-density 
(~ lO'' cm~^) outer parts of prestellar cores such as L1544 



( prutcher 1999| ; prutcher fc Troland 2000] ). Only relatively 
poor upper limits (< 100 /iG) exist for the field strength 
in the central parts of these cores (Levin et al. 2001). 
Albeit quite large the above value Bcrit ~ 60 /iG, which 
refers to the ^ 10^ cm~^ region, thus remains realistic. 
Furthermore, it should be noted that published ambipo- 
lar diffusion models include only a static magnetic field 
and do not take turbulent support into account. If a tur- 
bulent magnetic field is present, then a weaker static field 
may be sufficient to yield a decoupling density as high 
as ricrit ~ 10^ cm~^. According to this interpretation. 



Although rotating more slowly than the IRAM 04191 en- 
velope, some of the models are characterized by high ratios 
of centrifugal to gravitational acceleration at the supercritical 
radius, comparable to the 0.4 ± 0. 2 ratio o bserved here (cf. 
model 1 of BM94| and model 6 of |BM95b[ ). Qualitatively, a 
faster initial rotation rate is not expected t o change the evo- 
lution significantly compared to model 6 of BM95b (S. 
private communication). 



Basu, 
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Fig. 16. Specific angular momentum as a function of ra- 
dius for the Taurus 
al ~ 



ct 



NH3 dense cores of Goodman 
(19p3[) (open squares) and the rotating "envelopes" 
and "disks" observed by Ohashi and collaborators around 
Taurus Class I sources (filled circles and filled triangles, re- 
spectively) . The horizontal dashed line marks the cons tant 
specific angular momentum obtained by Ohashi (1999| ) for 
his sample of Class I sources; the solid line, correspond- 
ing to a constant angular velocity r2 = 1.8 km s~^ pc~^, 
shows the best-fit linear correlation determine d by us for 
the Taurus NH3 cores of Goodman ct al. (1993 ) . The loca- 
tions of the fastly rotating inner envelope (r ^ 3500 AU) 
and more slowly rotating outer envelope (r ^ 11000 AU) 
in IRAM 04191 are marked by open stars, based on the 
present study. (Adapted from Ohashi 1999| .) 



the observed radius of "decoupling" ~ 3500 AU would 
correspond to the cutoff wavelength for MH D waves, i.e., 
~ 6200 AUx (^) X ( 3>,iXm-3 )'' (cf- iMouschovias 



199l|)7|and the IRAM 04191 dense core would have ini- 
tially formed through the dissipation of MHD turbulence 



fusion. 

6.6. Evolution of angular momentum during protostar 
formation 

It is instructive to compare the rotational properties of 
the IRAM 04191 envelope with the characteristics of other 
prestellar and protostellar objects in Taurus. Our 2D ra- 
diative transfer modeling indicates a rotational velocity 
f rot = 0.20±0.04 km s~^ at r — 3500 AU, using an inclina- 
tion angle of z = 50° (cf. § 5.2). This corresponds to a local 



specific angular momentum j — 3.4 x 10 km s pc at 
Tin = 3500 AU (roughly equal to the mean half-power ra- 
dius of the N2H+ dense core: FWHM - 8800 AU x 5200 
AU) . Quite remarkably, the mean specific angular momen- 
tum J/M ~ 1.5 x 10""^ km pc we measure inside the 
rapidly rotating inner r < Tm envelope of IRAM 04191 

pc 



is very similar to the nearly constant ~ 10 km s ^ 



value found by Ohashi (1999) for a number of small-scale 
"envelopes" and "disks" around Taurus Class I sources. 



More precisely, in the diagram of specific angular momen- 
tum versus radius presented by phashi et al. (199"7b| ) (see 
their Fig. 6), IRAM 04191 hes close to the intersection 
between the "dense core" regime, where the angular ve- 
locity is approximately locked to a constant background 
value ~ 1 — 2 km pc~^ (presumably as a result of 
magnetic braking - see BM94| and § above) , and the 



"protostellar" regime, where the specific angular momen- 
tum is roughly constant in timeQ (cf. Fig. [l^ ). According 
to phashi et al. (19971 ), the transition between these two 
regimes at a radius ~ 5000 AU (i.e., ~ 0.03 pc) charac- 
terizes the size scale for dynamical collapse, inside which 
evolution proceeds with near conservation of angular mo- 
mentum. Interestingly, this size scale is comparable to the 
radius r^n ~ 3500 AU found here for the rapidly rotat- 
ing inner envelope of IRAM 04191. Our suggestion that 
the inner envelope is a magnetically supercritical core de- 
coupling from a subcritical environment (§ 6.5 above) is 
thus fully consistent with the finding and interpretation 
of [Ohashi et al. (1997b| ). 



7. Summary and conclusions 

We have carried out a detailed study of the structure and 
kinematics of the envelope surrounding the Class pro- 
tostar IRAM 04191 in Taurus. Our main results and con- 
clusions are as follows: 



1, 



(e.g.|Nakano 1998; Myers 1998) rather than ambipolar dif- 2. 



Extended, subsonic infall motions with Vi^f ^ 0.5 Os ~ 
0.1 km s^^, responsible for a marked 'blue infall asym- 
metry' in self-absorbed CS and II2CO lines, are present 
in the bulk of the envelope, up to at least ri^o ~ 
10000 — 12000 AU. The observations are also consis- 
tent with larger infall velocities scaling as finf oc r~^'^ 
in an inner region of radius ri « 2000 AU. The cor- 
responding mass infall rate is estimated to be Mjnf ^ 
2 - 3 X al/G ~ 3 X lO''^ Mq yr-\ 
The protostellar envelope is differentially rotating with 
an angular velocity profile Q oc 7--2-5±o.5 be^-^^een r-m « 
3500 AU and r^.o ^ 7000 AU. The rotation profile is 
shallower, albeit more poorly constrained, in the inner 

The angular velocity 



-0.9±0.4 



pc" 

'1 



at r ~ 3500 AU 



r < Tin region, i.e., oc r 
is estimated to be 51 ~ 12 km 

and only fl < 0.5 - 1 km s^^ pc~^ at r - 11000 AU 
The present value of the centrifugal radius is estimated 
to be less than 400 AU. 

The extended infall velocity prof ile is inconsistent with 
the in side-out collapse picture of ^hu, Adams fc Lizanc 



(1987 ) and only marginally consistent with isother- 
mal collapse models starting from marginally stable 
equilibrium Bonnor-Ebert spheres. The latter tend to 



^ Assuming that all Taurus sources follow a similar time evo- 
lution. Fig. hqmay be viewed as an evolutionary diagram where 
the radius plotted on the x-axis represents the "contraction 
state" of a core as a function of time. This is conceptually dif- 
ferent from a plot showing the spatial distribution of angular 
momentu m / velocity in a given core at a given time (such as 
Fig. Ic of Basu 1997 or Fig. [12|d of the present paper). 
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produce somewhat faster infall velocities than are ob- 
served. 

The contrast observed between the (steeply declining) 
rotation velocity profile and the (flat) infall velocity 
profile beyond ~ 3500 AU suggests that angu- 
lar momentum is not conserved in the outer envelope. 
This is difficult to account for in the context of non- 
magnetic collapse models. 

Based on a qualitative comparison with magnetic am- 
bipolar diffusion models of cloud collapse (e.g. BM94), 
we propose that the rapidly rotating inner envelope 
of IRAM 04191 corresponds to a magnetically super- 
critical core decoupling from an environment still sup- 
ported by magnetic fields and strongly affected by 
magnetic braking. In this view, the outer {r^ < r < 
i^ni,o) envelope represents a transition region between 
the forming protostar and the slowly rotating ambient 
cloud. Although published ambipolar diffusion models 
have difficulty explaining supercritical cores as small 
as Rcrit ~ 3500 AU, we speculate that more elabo- 
rate versions of these models, including the effects of 
MHD turbulence in the outer envelope, would be more 
satisfactory. 

Interestingly, the steepening of il{r) in IRAM 04191 
occurs at a radius comparable to the ^ 5000 AU scale 
inside which the specific angular momentum of Taurus 
dense cores appears to be conserved (cf. phashi et 
al. 1997b and Fig, p^ . Our results therefore support 
phashi et al. (1997b| )'s proposal that r - 5000 AU 



on Radiative Transfer in Molecular Lines held in Leiden in 
May 1999 (http://www.strw.leidenuniv.nl/~radtrans/). 
These tests correspond to a low-abundance and high- 
abundance HCO+ 12-level problem, respectively, in the 



represents the typical size scale for dynamical collapse 
in Taurus. More generally, we suggest that the rota- 
tion/infall properties observed here for IRAM 04191 
are representative of the physical conditions prevailing 
in isolated protostellar envelopes shortly (~ lO'* yr) 
after point mass formation. 
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Appendix: Characteristics of the MAPYSO 
code 

The numerical code we have used first calculates the 
non-LTE level populations with a ID (spherical) Monte- 
Carlo method (Bcrncs 1978, 1979). Radiative transfer 



along each line of sight and convolution with the antenna 
beam, approximated by a Ga ussian, are th en computed 
with the MAPYSO package ([Blinder 1997| ). The latter 
works in both ID and 2D geometry. 

We have tested the Monte Carlo code for two test prob- 
lems (1 and 2) available on the web page of the workshop 



context of the Shu (1977) spherical collapse model. The 
level populations and the excitation temperatures com- 
puted by our Monte-Carlo code without any reference field 
(see Bernes 197£; Pagani 1998) agree quite well with those 
calculated by the workshop participants. The only signif- 
icant difference occurs for the high-abundance case in the 
central region (inside ~ 1300 AU) where our HC0+(2- 
1), (3-2), and (4-3) excitation temperatures are lower by 
~ 20% compared to the main group results. 

The CS and C'^^'S Monte-Carlo calculations reported 
in § ^ and § || used 9 levels and 27 concentric shells. This 
number of levels should be sufficient as the 9"" level is 
85 K above the ground level while the kinetic tempera- 
ture in the envelope does not exceed 20 K (see § LI). We 
used the CS collision rates computed in the 20-300 K range 
by Turner et al. (1992) and extrapolated these to 5-300 K 
with polynomials ( phoi et al. 1995 and N. Evans, private 
communication). Each simulation was performed with- 
out any reference field and resulted from two successive 
Monte-Carlo runs. The first run started from LTE, used 
packets of 1000 model photons, and computed 100 itera- 
tions, reinitializing the counters after each iteration (see 



Bernes 1979). It converged rapidly but still suffered from 



a high level of statistical noise. The second run improved 
the convergence and reduced the noise level by starting 
from the output of the first run, computing 40 iterations 
with packets of 40000 model photons, and reinitializing 
the counters after each group of five iterations. We checked 
that this number of iterations was large enough to reach 
convergence on the populations of the first five levels with 
an accuracy better than a few percents. 



References 

Alves, J.F., Lada, C.J., & Lada, E.A. 2001, Nature, 409, 
159 

Andre, P., Motte, F., & Bacmann, A. 1999, ApJ, 513, L57 
(AMB99) 

Andre, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 
406, 122 

Andre, P., Ward-Thompson, D., & Barsony, M. 2000, in 
Protostars and Planets IV, ed. V. Mannings, A. P. Boss, 
& S.S. Russcl (Univ. of Arizona Press, Tucson), p. 59 

Andre, P., Ward- Thompson, D., & Motte, F. 1996, A&A, 
314, 625 

Arquilla, R., & Goldsmith, P.F. 1986, ApJ, 303, 356 
Artymowicz, P., & Lubow, S.H.1994, ApJ, 421, 651 
Bachiller, R., & Prez Guttirrez, M. 1997, ApJ, 487, L93 
Bacmann, A., Andre, P., Puget, J.-L., Abergel, A., 

Bontemps, S., & Ward-Thompson, D. 2000, A&A, 361, 

555 

Basu, S. 1997, ApJ, 485, 240 
Basu, S. 1998, ApJ, 509, 229 



22 



Belloche et al.: Velocity structure of the IRAM 04191 envelope 



Basu, S., & Mouschovias, T.C. 1994, ApJ, 432, 720 
(BM94) 

Basu, S., & Mouschovias, T.C. 1995a, ApJ, 452, 386 
(BM95a) 

Basu, S., & Mouschovias, T.C. 1995b, ApJ, 453, 271 

(BM95b) 

Beckwith, S.V.W., Sargent, A.I., Chini, R.S., & Guesten, 

R. 1990, AJ, 99, 924 
Belloche, A., Andre, P., & Motte, F. 2001, in From 

Darkness to Light, Eds. T. Montmerle & P. Andre, ASP 

Conf. Ser., 243, p. 313 
Benson, P.J., & Myers, P.C. 1989, ApJSS, 71, 89 
Bergin, E.A., & Langer, W.D. 1997, ApJ, 486, 316 
Bernes, C. 1978, Stockholm Observatory Report no. 15 + 

erratum 
Bernes, C. 1979, A&A, 73, 67 
Blinder, S. 1997, Thesis, Universite de Bordeaux I 
Bonnor, W.B. 1956, MNRAS, 116, 351 
Broguiere, D., Neri, R., Sievers, A., & Wiesemeyer, H. 

2002, NIC Bolometer Users Guide, V 1.4-11 
Buisson, G., Desbats, L., Duvert, G., Forveille, T., Gras, 

R., Guilloteau, S., Lucas, R., & Valiron, P. 2002, CLASS 

manual, Version 4.0 (IRAM) 
Burkert, A., & Bodenheimer, P. 2000, ApJ, 543, 822 
Cabrit, S., & Bertout, C. 1990, 348, 530 
Caselh, P., Myers, P.C, & Thaddeus, P. 1995, ApJ, 455, 

L77 

Caselli, P., Walmslcy, CM., Terzevia, R., & Herbst, E. 

1998, ApJ, 499, 234 
Caselli, P., Walmslcy, CM., Zucconi, A., Tafalla, M., 

Dore, L., & Myers, P.C. 2001a, ApJ, 565, 331 
Caselh, P., Walmslcy, CM., Zucconi, A., Tafalla, M., 

Dore, L., & Myers, P.C. 2001b, ApJ, 565, 344 
Ceccarelh, C, HoUcnbach, D.J., & Tielens, A.G.G.M. 

1996, ApJ, 471, 400 
Chevalier, R.A. 1983, ApJ, 268, 753 

Choi, M., Evans, N.J. 11, Grcgcrscn, E., & Wang, Y. 1995, 

ApJ, 448, 742 
Ciolck, G.E., & Basil, S. 2000, ApJ, 529, 925 
Ciolek, G.E., & Kmigl, A. 1998, ApJ, 504, 257 
Ciolek, G.E., & Mouschovias, T.C. 1994, ApJ, 425, 142 
Crutcher, R.M. 1999, ApJ, 520, 706 
Crutcher, R.M., & Troland, T.H. 2000, ApJ, 537, L139 
Doty, S.D., & Ncufcld, D.A. 1997, ApJ, 489, 122 
Evans, N.J. II 1999, ARAA, 37, 311 
Evans, N.J. II, Rawlings, J.M.C, Shirley, Y.L., & Mundy, 

L.G. 2001, ApJ, 557, 193 
Foster, P.N., & Chevalier, R.A. 1993, ApJ, 416, 303 
Galh, D., & Shu, F.H. 1993, ApJ, 417, 220 
Goodman, A. A., Benson, P.J., Fuller, G.A., & Myers, P.C. 

1993, ApJ, 406, 528 
Gottlieb, C.A., Myers, P.C, & Thaddeus, P. 2002, in 

preparation 

Guilloteau, S., Lucas, R., & Bouyoucef, K. 2002, Mapping 
Cookbook IRAM Plateau de Bure Interferometer, 
Version 1.3 (IRAM) 

Hartmann, L. 1998, Accretion Processes in Star 
Formation, Cambridge Astrophysics Series, Vol. 32 



(Cambridge University Press) 
Hennebelle, P., Whitworth, A.P., Gladwin, P.P., & Andr, 

P. 2002, MNRAS, in press 
Henriksen, R.N., Andre, P., & Bontemps, S. 1997, A&A, 

323, 549 (HAB97) 
Hogerhcijde, M.R., & SandeU, G. 2000, ApJ, 534, 880 
Jensen, E.L.N., Mathieu, R.D., & Fuller, G.A. 1996, ApJ, 

458, 312 

Jessop, N.E., & Ward-Thompson, D. 2000, MNRAS, 311, 
63 

Kenyon, S.J., Calvet, N., & Hartmann, L.W. 1993, ApJ, 

414, 676 

Klapper, G., Lewen, F., Gendriesch, R., Belov, S.P., & 
Winnewisser, G., 2001, Z. Naturforsch., 56a, 329 

Kramer, C, Stutzki, J., Rohrig, R., & Corneliussen, U. 
1998, A&A, 329, 249 

Krasnopolsky, R., & Konigl, A. 2002, ApJ, submitted 

Larson, R.B., 1969, MNRAS, 145, 271 

Lee, C.W., Myers, P.C, & Tafalla, M. 2001, ApJS, 136, 
703 

Levin, S.M., Langer, W.D., Velusamy, T., Kuiper, T.B.H., 

& Crutcher, R.M. 2001, ApJ, 555, 850 
Li, Z.-Y. 1998, ApJ, 493, 230 
Li, Z.-Y., & Shu, F.H. 1996, ApJ, 472, 211 
Li, Z.-Y., & Shu, F.H. 1997, ApJ, 475, 237 
Lucas, R. 1999, CLIC manual. Version 4.1 (IRAM) 
Lucas, R. 2001, in Proceedings from IRAM Interferometry 

School 2, Ed. A. Dutrey (IRAM), p 175 
Masunaga, H., & Inutsuka, S.I. 2000, ApJ, 531, 350 
Masunaga, H., Miyama, S.M., & Inutsuka, S. 1998, ApJ, 

495, 346 

Motte, F., & Andr, P. 2001, A&A, 365, 440 (MAOl) 
Motte, F., Andre, P., & Neri, R. 1998, A&A, 336, 150 

(MAN98) 

Mouschovias, T.C. 1991, ApJ, 373, 169 
Myers, P.C. 1998, ApJ, 496, L109 

Myers, P.C. 1999, in The Physics and Chemistry of the 
Interstellar Medium, Proceedings of the 3rd Cologne- 
Zermatt Symposium, Eds : V. Ossenkopf. J. Stutzki, 
and G. Winnewisser (GCA-Verlag Herdecke) 

Myers, P.C, Bachillcr, R., CaseUi, P., Fuller, G.A., 
Mardoncs, D., Tafalla, M., & Wilner, D.J. 1995, ApJ, 
449, L65 

Myers, P.C, Evans, N.J. II, & Ohashi, N. 2000, in 
Protostars and planets IV, Ed. V. Mannings, A. Boss, 
& S. Russel (Tucson: Univ. Arizona Press), p. 217 

Nakano, T. 1998, ApJ, 494, 587 

Ohashi, N. 1999, in Star Formation 1999, Ed : T. 

Nakamoto (Nobeyama Radio Observatory), pl29 
Ohashi, N., Hayashi, M., Ho, P.T.P., & Momose, M. 1997a, 

ApJ, 475, 211 

Ohashi, N., Hayashi, M., Ho, P.T.P., Momose, M., 
Tamura, M., Hirano, N., & Sargent, A.I. 1997b, ApJ, 
488, 317 

Onishi, T., Mizuno, A., & Fukui, Y. 1999, PASJ, 51, 257 
Pagani, L. 1998, A&A, 333, 269 
Penston, M.V., 1969, MNRAS, 144, 425 
Saigo, K., & Hanawa, T. 1998, ApJ, 493, 342 



Belloche et al.: Velocity structure of the IRAM 04191 envelope 



23 



Safier, P.N., McKee, C.F., & Stabler, S.W. 1997, ApJ, 485, 
660 

Shu, F. 1977, ApJ, 214, 488 

Shu, F.H., Adams, F.C., & Lizano, S. 1987, ARA&A 25, 
23 

Spitzer, L., Jr. 1978, Physical Processes in the Interstellar 

Medium (New York: Wiley) 
Stabler, S.W., Korycansky, D.G., Brothers, M.J., & 

Touma, J. 1994, ApJ, 431 341 
Stutzki, J. , & Glisten, R. 1990, 356, 513 
Tafalla, M., Mardones, D., Myers, P.C., Caselli, P., 

Bachiller, R., & Benson, P.J. 1998, ApJ, 504, 900 
Tafalla, M., Myers, P.C., Caselli, P., Walmsley, CM ., & 

Comito, C. 2002, ApJ, in press ( |astro-ph/0112487D 
Tamura, M., Gatley, I., Waller, W., & Werner, M.W. 1991, 

ApJ, 374, L25 

Terebey, S., Shu, F.H., & Cassen, P. 1984, ApJ, 286, 529 
(TSC84) 

Tomisaka, K. 1996, PASJ, 48, L97 
Tomisaka, K. 2000, ApJ, 528, L41 

Turner, B.E., Chan, K.W., Green, S., & Lubowich, D.A. 

1992, ApJ, 399, 114 

Ungerecbts H., et al. 2000, in Imaging at Radio through 
Submilbmeter Wavelengths, ASP Conf. Proc, Vol. 217, 
Eds : J.G. Mangum & S.J.E. Radford (Astronomical 
Society of the Pacific), pl90 

Vrtilek, J.M., Gottlieb, C.A., & Thaddeus, P. 1987, ApJ, 
314, 716 

Walker, C.K., Narayanan G., & Boss, A.P. 1994, ApJ, 431, 
767 

Ward-Thompson, D., & Buckley, H.D. 2001, MNRAS, 327, 
955 

Ward-Thompson, D., Motte, F., & Andre, P. 1999, 

MNRAS, 305, 143 
Ward-Thompson, D., Scott, P.F., Hills, R.E., & Andre, P. 

1994, MNRAS, 268, 276 
Whitworth, A., & Summers, D. 1985, MNRAS, 214, 1 
Wild, W. 1999, The 30m Manual, Pico Veleta, Spain, 

Version 2.0 (IRAM) 
Williams, J.P., & Myers, P.C. 1999, ApJ, 518, L37 
Wilner, D.J., Myers, P.C, Mardones, D., & Tafalla, M. 

2000, ApJ, 544. L73 
Wilson, T.L., & Rood, R.T. 1994, ARA&A, 32, 191 
Wootten, A., Wiseman, J., & Fuller, G. 2001, AAS 

Meeting 199, #60.02 
Zhou, S. 1995, ApJ, 442, 685 

Zhou, S., Evans, N.J. II, Kompe, C, & Walmsley, CM. 

1993, ApJ, 404, 232 

Zhou, S., Evans, N.J. II, & Wang, Y. 1996, ApJ, 466, 296 
Zucconi, A., Walmsley, CM., & Galh, D., 2001, A&A, 376, 
650 



